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Preface 


During  the  past  year  the  NRL  Advanced  Simulation  Concepts  Theory 
Program  has  focused  on  extending  and  improving  the  non-LTE  equation  of 
state  data  base,  investigating  the  dynamics  of  imploding  Z-pinches,  with 
and  without  a  central  plasma  core,  predicting  the  radiation  from  ion  beam 
heated  targets,  developing  a  simple  K-shell  scaling  law  for  conventional 
PRS  loads,  and  the  preparation  of  a  pulse  power  primer  for  the  beginner. 
The  non-LTE  equation  of  state  work  was  devoted  to  studying  the  radiation 
enrichment  of  gas  puff  mixtures,  i.e.,  determine  whether  it  is  possible  to 
optimize  the  power  radiated  by  a  plasma  comprised  of  two  elements.  Two 
cases  were  investigated:  the  optically  thin  and  optically  thick  plasmas, 
respectively.  For  the  optically  thin  investigation  the  plasma  was  composed 
of  two  elements  of  widely  separated  atomic  number.  The  results  suggest 
that  the  efficient  radiation  from  lower  shells  of  higher-Z  elements  require 
more  powerful  generators  to  attain  the  same  conditions  as  higher  3hells  of 
lower-Z  loads.  A  not  too  surprising  conclusion.  However,  it  was  3hown 
that  "spectral  fill-in"  can  be  obtained  by  mixtures  using  appropriate 
concentrations  of  complementary  constituents,  resulting  in  smoother  photon 
distributions  in  the  line  spectrum.  For  the  optically  thick  scenario  the 
plasma  was  composed  of  two  elements  with  adjacent  atomic  number.  Under 
these  conditions  the  results  show  convincingly  that  the  K-shell  emissivity 
of  the  plasma  can  be  increased  signif iciantly  over  the  conventional  single¬ 
element  load.  The  mechanism  that  increases  the  emissivity  is  that  the 
opacity  of  each  element  is  less,  allowing  more  photons  out  -  and  twice  as 
many  lines  are  present  in  the  K-shell  region. 

Numerical  simulations  were  performed  to  provide  guidance  for  the  gas 
puff  experiments  conducted  on  the  GAMBLE  II  facility  and  to  determine  the 
stability  of  plasma  shells  during  the  runin  acceleration  phase  with  Double- 
Eagle  class  generators.  The  GAMBLE  II  experiments  were  modeled  with  the 
1-D  radiation  hydrodynamics  simulation  code,  PRISM.  The  numerical 
simulations  compared  and  contrasted  the  differences  in  code  results  with 
and  without  radiation  transport  and  detailed  atomic  physics  modeling.  The 
code  comparisons  yield  larger  differences  than  those  obtained  from  a 
comparison  of  the  full  model  and  the  experimental  results.  The  non-LTE 
numerical  simulations  were  in  particularly  good  agreement  with  the 
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experimental  results.  On  the  question  of  plasma  stability  during  the  run- 
in  acceleration  phase  two  issues  were  considered:  first,  the  amount  of 
shell  mass  participating  in  the  unstable  configuration  and  second,  the 
magnitude  of  the  growth  rate  of  the  instability  and  how  it  compares  with 
the  implosion  time.  It  was  determined  that  about  60?  of  the  initial  mass 
participates  in  the  implosion  to  completion  and  although  the  back  surface 
starts  breaking-up,  most  of  the  mass  remains  stable  and  can  become  the 
source  of  intense  radiation  i.e.,  degradation  of  K-line  radiation  is  only 
minimally  influenced  by  the  instability. 

Another  area  of  considerable  interest  is  the  examination  of  stagnating 
an  imploding  gas  puff  plasma  onto  a  low  density  plasma  for  both  radiation 
enhancement  or  for  use  as  a  fast  switch.  Application  of  plasma  stagnation 
as  a  fast  switch  is  currently  being  examined.  For  the  purposes  of  this 
report  the  process  of  plasma  stagnation  is  investigated  as  a  source  of 
radiation  enhancement.  The  numerical  simulations  deal  with  the  implosion 
of  a  neon  gas  puff  onto  a  central  neon  core.  The  densities  and 
temperatures  at  stagnation  over  interesting  volumes  and  timescales  are 
sufficient  to  suggest  the  efficient  conversion  of  kinetic  to  thermal  energy 
which  is  subsequently  converted  to  radiation.  About  two-thirds  of  the 
total  plasma  energy  was  converted  to  K-  and  L-shell  radiation  in  100 
nanoseconds. 

Velocity  scaling  of  a  radiatively  dominated  Z-pinch  was  investigated 
using  an  electrodiffusive  model.  Preliminary  results  show  that  the  effect 
of  the  early  radiation  on  the  scaling  error  is  most  acute  for  the  low 
current  cases.  However,  the  general  trend  is  to  reduce  the  pinch  velocity 
and  narrow  the  load  width.  A  more  extensive  study  is  in  progress. 

Finally,  an  improved  model  for  evaluating  proton  stopping  power  was 
developed  and  incorporated  into  the  hydrodynamic  response  simulation  code 
in  place  of  the  Local  Oscillator  Model  (LOM).  The  new  version  incorporates 
the  results  of  detailed  theoretical  calculations  of  proton  stopping  cross 
sections  by  A1  ions  from  a  Generalized  Oscillator  Strength  treatment  of  the 
projectile  target  interaction.  For  neutral  aluminum  the  agreement  with  the 
LOM  result  is  good,  but  predicts  smaller  stopping  power  than  LOM  results 
for  ionized  A1 ,  particularly  at  lower  proton  energies  and  higher  charge 
states.  These  lower  stopping  powers  result  in  deeper  penetration  into  the 
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target  by  the  proton  beam,  especially  at  early  times.  This  causes  the  rear 
side  to  be  hotter  and  less  dense  than  if  LOM  stopping  powers  are  used.  At 
times  of  the  order  of  a  nonosecond,  the  rear  side  plasma  emission  is 
greater  than  the  front  side;  when  the  LOM  is  used,  the  reverse  situation 
occurs.  Hence,  the  use  of  different  stopping  powers  can  result  in  marked 
differences  in  plasma  emission  spectra. 

In  addition  to  the  work  presented  in  this  report  there  has  been  a 
number  of  reports  on  related  activities.  They  include: 

"Imploding  Plasma  Radiation  Sources:  Basic  Concepts." 

NRL  Memo  Rept.  5392  (July,  1984) 

"K-Shell  Yield  Scaling  Law  for  Conventional  PRS  Loads." 

NRL  Memo  Rept.  5406  (Aug.  1984) 

"Hydrodynamic  Response  of  a  Radiating  A1  Plasma  Excited 
by  a  Proton  Beam."  NRL  Memo  Rept.  5482  (Dec.  1984). 

Copies  of  these  reports  are  available  on  request. 
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A.  Radiative  Properties  of  Puffed-Gas  Mixtures:  Optically  Thin  Case 


I.  INTRODUCTION 

It  has  been  a  major  goal  of  the  plasma  physics  community  in  recent 

years  to  produce  a  high-fluence  source  of  soft  x-rays  (hv>  1.0  keV)  in  the 

laboratory.  The  applications  of  such  a  source  are  numerous:  x-ray 

lithography,  x-ray  microscopy,  extended  x-ray  absorption  fine  structure 

(EXAFS)  studies,  photon-pumped  short-wavelength  lasers,  and  x-ray 
vulnerability  studies,  to  name  a  few  examples.  In  the  effort  to  construct 
such  a  source,  several  diverse  experimental  configurations  have  been 
successfully  demonstrated,  such  as  laser-target  and  ion  beam-target 
interaction,  synchrotons ,  imploding  foil  liners,  imploding  wire  arrays ,  and 
gas-puff  Z-pinches.  The  last  of  these  techniques1-10  rely  on  high-voltage, 
pulse-power  discharges  through  a  plasma  load,  transforming  the  low- 
temperature  material  into  an  imploding  annular  plasma,  accelerated  radially 
by  the  current-dr i ven ,  azimuthal  magnetic  field  to  velocities  near  107 
cm/sec.  The  on-axis  assembly  of  the  imploding  plasma  results  in 
thermalization  of  the  kinetic  energy  driving  it  and  the  subsequent  emission 
of  copious  radiation  in  the  soft  x-ray  band.  Due  to  the  electrical 
character istics  of  these  machines  (e.g.,  attainable  peak  voltage  and 
current,  impedance,  inductance,  and  pulse  width),  much  time  and  effort  has 
been  invested  in  matching  the  appropriate  load  (be  it  a  foil,  wires  or  gas- 
puff)  to  the  generator  in  order  to  optimize  x-ray  emission  in  the  desired 
wavelength  regions.  One  proposed  method  of  enhancing  the  x-ray  yield  (as 
well  as  broadening  the  spectral  region  where  emission  occurs)  is  to  replace 
single-material  loads  with  mixtures.  Although  little  work  has  been  devoted 
to  study  this  proposed  technique,  '  it  has  been  suggested  that  by  combining 
materials  of  low  and  high  atomic  number,  an  increase  in  either  the  plasma 
temperature  or  density  (or  both)  during  the  x-ray  emission  phase  would  be 
realized,  a3  compared  to  that  of  purely  the  higher-Z  material  implosion. 

p  a 

Since  the  desired  soft  x-ray  emission  scales  nearly  as  Nj  and  up  to  Tg  in 
certain  regimes  (Nj  and  Te  being  plasma  ion  density  and  electron 
temperature),  it  is  hoped  that  the  increase  in  these  parameters  would  yield 
an  increment  in  radiation  which  would  exceed  that  which  was  lost  by  the  use 
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f  i  I  ewer  concentration  of  the  hi gher-Z  material,  thereby  enhancing  the 
-’aoiative  output. 

In  this  study,  a  theoretical  investigation  of  the  radiative  proper¬ 
ties  of  multi-materi  al  gas-puff  Z-pinch  plasmas  is  described.  The  primary 
tool  used  in  this  study  is  a  numerical  model  of  plasma  ionization 
/radiation  with  the  newly-developed  capability  of  treating  multi-material 
plasmas  simultaneously  and  self-consistently .  The  model,  which  is 
discussed  below  in  more  detail,  is  based  on  a  solution  of  a  set  of  atomic 
rate  equations  characterizing  a  plasma  in  coll isional-radi ati ve  equilibrium 
(CRE).  °  Parameter  studies  are  presented  for  mixtures  consisting  of 
various  concentrations  of  helium/argon,  hel ium/ krypton ,  neon/argon,  and 
argon/ krypton;  the  radiative  yields  are  reported  as  a  function  of  plasma 
temperature  '~y  maintaining  a  constant  mass  density  of  the  mixture.  In 
order  to  separate  the  purely  collisional  effects  on  the  radiation 
enhancement  from  photon  pumping  processes,  only  optically  thin  plasmas  are 
considered  here.  The  cooperative  effects  of  collisions,  opacity  and 
radiation  transport  will  be  studied  in  a  paper  to  appear  later. 

I!.  MODEL,  DESCRIPTION 

a.  loni zation/Radiation  Model 

The  basic  ionization/radiation  model  has  been  discussed  and 

1  Q 

documented  in  several  previous  publications  7  and  only  a  brief  description 
will  be  given  here.  A  set  of  atomic  rate  equations  of  the  form, 

dN. 

~  «  £  W..N,  -  Z  W. .N.  (1) 

dt  J  ji  j  J  ij  i 

are  solved  for  the  ground  and  excited  state  population  densities  N^,  where 
W  j  j  is  a  collisional  or  radiative  transition  rate  from  state  i  to  state 
j.  The  processes  comprising  Wy  used  in  this  investigation  include: 
collisional  ionization,  collisional,  radiative,  and  dielectronic 
recombination,  collisional  excitation  and  deexcitation,  and  spontaneous 
radiative  decay.  As  mentioned  above,  this  study  assumes  the  plasmas  to  be 
optically  thin,  hence,  photoioni zat ion ,  photoexcitation  and  stimulated 
emission  processes  are  neglected.  The  rate  coefficients  used  to  describe 
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the  plasma  atomic  transitions  are  documented  in  the  next  section.  Note 
that  this  is  an  investigation  of  the  radiative  properties  of  mixtures  in 
equilibrium  at  some  nominal  electron  temperature  and  ion  density.  The 
combined  effects  of  the  radiation  and  the  magnetohydrodynamics  (MHD)  of  the 
current-dri ven  plasma  mixture  are  a  more  complicated  issue  and  will  be 
addressed  at  a  later  time. 

In  order  to  treat  the  dual-component  plasma  mixture,  the  set  of  atomic 
rate  equations  for  each  material  must  be  solved  simultaneously,  since  the 
free  electrons  due  to  ionization  of  both  materials  are  shared  in  the  sense 
that  any  electron  can  collide  with  any  ion.  Thus,  the  electron  density  for 
a  given  temperature  and  fixed  total  mass  density  will  lie  between  that  of 
the  pure  lower-Z  plasma  and  the  pure  higher-Z  plasma,  and  is  given  by, 

N-EZ.N.  +  E  z!  1  (2) 

e  ,xl  .11 

l  l 

where  Z;  and  Z!  are  the  charges  of  state  i  in  material  1  and  material  2  (of 
a  l 

course,  the  level  structure  of  each  material  employed  in  the  model  may  be 
different).  Since  the  electron  density  occurs  in  Wy  (via  both  2-body  and 
3-body  processes),  the  set  of  equations  is  highly  non-linear  and  is  solved 
numerically  using  an  iterative  predictor-corrector  scheme. 

.  The  plasma  radiative  yield  is  comprised  of  emission  via  three  atomic 
mechanisms:  (i)  line  radiation  from  bound-bound  transitions,  (ii) 

recombination  radiation  from  bound-free  transitions,  and  (iii) 
bremsstrahlung  emission  via  free-free  interactions.  At  densities  and 
temperatures  typical  of  most  plasma  implosions  on  the  largest  pulse  power 
generators  presently  configured  for  gas  puffs  ("BLACKJACK  V"'  at  Maxwell 
Laborator ies ,  "PITHON"  at  Physics  International,  "SHIVA"  at  the  Air  Force 
Weapons  Laboratory,  "PROTO  II"  at  Sandia  Laborator ies ,  and  "GAMBLE  II"  at 
the  Naval  Research  Laboratory),  about  0. 5~2.0  keV  at  “lO1^  ions/cm^,  the 
largest  fraction  of  the  total  radiation  is  emitted  as  lines.  ’  For  this 
reason  it  is  necessary  to  take  into  account  all  of  the  strong  lines  emitted 
in  this  temperature  range  to  accurately  determine  the  yield.  The  level 
structure  employed  to  accomplish  this  is  described  below. 
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o.  Atomic  Data 


A  critical  factor  affecting  the  integrity  of  a  radiation 
calculation  is  both  the  accuracy  and  completeness  of  the  data  base  used  in 
the  ionization/radiation  calculation.  This  includes  the  level  energies, 
rate  coefficients  and  the  number  of  excited  levels  of  each  ion  that  are 
retained  in  the  model.  This  investigation  included  studies  with  the  first 
four  noble  gases:  helium,  neon,  argon,  and  krypton.  Since  they  vary  in 
atomic  number  from  2  to  36,  and  the  temperature  range  of  interest  was  held 
relatively  constant,  different  level  structure  was  required  for  each 
material.  Helium  was  assumed  to  be  fully  ionized  at  temperatures  above  100 
e'7  and,  hence,  only  the  bare  nucleus  and  two  free  electrons  are  free  to 
interact.  No  detailed  level  structure  was  therefore  needed  for  helium. 
Neon  will  most  likely  be  stripped  of  its  1-shell  electrons  above  100  eV, 
hence,  excited  level  structure  was  retained  for  the  lithiumlike, 
helium! ike,  and  hydrogenlike  ions.  Argon  will  have  lost  its  M-shell 
electrons  at  100  eV,1M  so  structure  was  included  for  Ar  IX-Ar  XVIII 
(neonlike  through  hydrogenlike).  Krypton,  on  the  other  hand  will  not  be 
stripped  down  to  the  L-shell  electrons  until  about  1.0  keV.^0  Since  at 
present  we  have  modeled  the  excited  level  structure  in  krypton  only  down  to 
Kr  XXVI  (sodiumlike  ion)  we  cannot  conduct  meaningful  calculations  of  total 
radiative  yield  below  temperatures  of  500  eV.  However,  for  our  purposes, 
we  are  predominantly  interested  in  photon  emission  above  a  kilovolt;  since 
the  M-shell  ions  emit  lines  below  about  900  eV  and  L-shell  ions  emit  most 
lines  (An21)  above  1.6  keV,  all  of  the  relevant  line  emission  is  included 
in  the  model.  A  list  of  the  energy  levels  included  in  the  model  is  given  in 
Table  I;  the  same  structure  was  used  for  both  argon  and  krypton,  while  only 
the  1-,  2-,  and  3-electron  systems  apply  to  neon.  All  ground  states  were 
included  for  each  of  the  three  materials. 

The  level  energies  and  radiative  decay  coefficients  for  the  level 
structure  in  Table  I  were  calculated  using  the  Hartree-Slater  atomic 

op 

structure  model  written  by  R.  Cowan.  Since  this  code  provides 

information  for  individual  j-levels,  significant  averaging  over  j  (and 
oftentimes  2. )  was  necessary  to  obtain  energies  and  Einstein-A  coefficients 
for  the  structure  given  in  Table  I.  Electron  impact  excitation  is 
calculated  by  a  number  of  techniques  and  a  short  discussion  of  these 


follows.  In  the  case  of  krypton,  a  complete  set  of  collision  strengths 
calculated  oy  Bhatia  and  Feldman--1  in  the  distorted  wave  approximation  were 
used  to  generate  the  rate  coefficients  for  the  An=0  transitions  in  L-shell 
ions.  The  rate  coefficients  for  K-shell  and  lithiumlike  transitions  in 
argon  were  calculated  using  the  method  of  distorted  waves  described  by 

oy 

Blaha,-  while  those  for  neon  and  krypton  were  calculated  using  the 
Coulomb-Borr.  scheme  in  intermediate  coupling  described  by  Sampson  and  co- 
workers."  The  screening  constants  (both  direct  and  exchange 

contributions)  used  in  this  calculation  were  obtained  via  direct 
comparisons  with  the  argon  distorted  wave  results.  The  transitions 

coupling  the  n  =  1  and  n=2  levels  in  the  heliumlike  ions  of  neon  and  argon 
were  taken  from  the  results  of  Pradhan  et  al .  Q  which  include  the  effects 
of  autoionizing  resonances.  The  rates  for  sodiumlike  krypton  were  taken 

O  'T 

from  the  distorted  wave  calculations  of  Blaha  and  Davis,  '  while  those  for 

pQ 

neor.like  krypton  have  been  calculated  by  Reed  and  Hazi.  Rate 

coefficients  for  carbonlike,  boronlike,  and  berylliumlike  argon  were  also 

calculated  using  31aha's  distorted  wave  model.  All  other  rate  coefficients 

not  obtainable  by  more  sophisticated  means  were  calculated  using  a  code 

pq 

developed  by  Davis  7  based  on  the  semi-classical  impact  parameter 
technique;  the  oscillator  strengths  calculated  using  the  atomic  structure 
code  were  used  as  the  input,  and  bound-bound  Gaunt  factors  were  calculated 
directly.  The  rate  coefficients  for  the  remaining  collisional  processes 
are  identical  to  those  used  in  previous  investigations  and  the  details  of 
the  calculations!  methods  can  be  found  by  consulting  Ref.  18. 

III.  RESULTS 

a.  Density  Scaling 

As  was  stated  earlier,  no  attempt  was  made  to  account  for  the 
dynamics  of  the  imploding  plasma  in  tnis  study,  hence,  the  radiation/ 
ionization  calculation  took  electron  temperature  and  the  density  of  each 
constituent  as  input  variables.  In  this  sense,  the  predicted  radiative 
yields  in  a  given  spectral  range  are  representati ve  of  only  the  hottest, 
most  dense  part  of  the  pinched  plasma,  along  the  axis  of  the  initial  gas 
annulus.  3revious  MHD  simulations^  indicate  that  surrounding  the  hotter 
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plasma  are  colder,  less  dense  plasma  regions  character i zea  by  large 
temperature  gradients.  Since  the  ions  found  in  these  cooler  regions  emit 
radiation  at  frequencies  below  the  range  of  interest  here,  it  will  suffice 
to  neglect  these  regions  and  treat  only  the  hotter  plasma.  The  peak 
electron  temperatures  attained  in  these  plasmas  usually  vary  according  to 
the  size  of  tne  pulse  power  generator  and  the  initial  gas  pressure. 
However,  the  '-ange  of  interest  is  typically  between  300  eV  and  2.5  keV, 
therefore,  we  limited  our  parameter  study  to  this  region. 

? cur  different  gas  mixtures  were  studied  for  radiation  enhancement: 
neon/ argon,  argon/ krypton,  helium/ argon ,  and  heli urn/ krypton .  The  spectral 
ranges  of  interest  were  above  900  eV  for  neon  (K-snell  emission),  above  3.0 
keV  for  argon  (X- shell  emission),  and  above  1.5  keV  for  krypton  (L-shell 
emission).  Helium  emits  only  icw-fluenee  fcremsstrahlung  in  these  spectral 
ranges  and  was  therefore  neglected  in  the  yield  analysis.  We  conducted  the 
parameter  study  by  varying  the  concentration  of  the  gases  in  each  mixture 
for  a  number-  of  relevant  temperatures.  The  present  level  of  understanding 
of  gas- puff  implosions  dictates  that  -when  varying  the  concentrations  of  the 
mixture  gases,  one  should  hold  the  total  mass  density  of  the  mixture 
constant  in  order  to  make  fair  comparisons  of  the  yield.  In  this  way,  as 
the  concentrations  vary,  the  same  mass  will  be  accelerated,  presumably 
attaining  the  same  temperature  at  assembly.  Thus  the  ion  density  of 
constituent  1  is  given  by 


m 


an:  that  of  constituent  2  by 


m 


<«) 


where  m7  is  the  mass  density  in  amu/cm^,  a,  and  A2  are  the  atomic  weights 
of  the  constituents  and  a  is  the  ion  number  concentr ation  of  constituent 
.  d  =  'I ,  ','J.  *  Ml).  The  mass  density  was  held  fixed  at  2x10^®  amu/cm^ 

1  Q  V 

wn.cn  is  commensurate  with  a  density  of  10  7  ions /cm-*  for  a  pure  neon 
clasma.  This  corresponds  approximately  to  the  densities  determined  from 
spectral  analysis  of  K-shell  emission  on  the  Blackjack  V  machine,  although 
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Total  radiation  versus  argon  concentrat i on  for  neon/ argon 
mixtures;  solid  curves  are  argon  *  neon,  broken  curves  are  the 


a rgon  alone. 


log-jQ  YIELD  (W/crrr) 


|J I  I  lip  III  I  H  IJI  ■  W  M  j»  in  lin 


Fig.  3.  Line  radiation  yield  (hv  >  1.5  keV)  versus  krypton  concentration 
for  ar gon/ kr ypton  mixtures;  solid  curves  are  argon  +  krypton, 
broken  curves  are  the  krypton  contribution. 
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%  CONCENTRATION  ARGON 


Line  radiation  yield  (hv  >0.9  keV)  versus  argon  concentration 
for  neon/argon  mixtures;  solid  curves  are  argon  +  neon,  broken 
curves  are  the  argon  contribution. 
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Table  I  (Cont’d)  —  Atomic  level  structure  used  in  the  ionization/radiation  model  for  neon, 
argon,  and  krypton  (only  the  krypton  model  has  Na-like  excited  levels). 

Li -like 


ls22s 

2 

2p 

6 

3s 

2 

3p 

6 

3d 

10 

42 

32 

52 

50 

He- like 


Is2  ( 

1S)  1 

ls2s 

(3S) 

3 

ls2s 

(LS) 

1 

ls2p 

(3P) 

9 

ls2p 

(XP) 

3 

Is  32 

(3l) 

27 

ls32 

(ll) 

9 

Is  42 

64 

ls52 

100 

H-like 


Table  I  (Cont’d)  —  Atomic  level  structure  used  in  the  ionization/radiation  model  for  neon, 
argon,  and  krypton  (only  the  krypton  model  has  Na-like  excited  levels). 


N-like 

B-like 

ls22s22p3  (4S3/2) 

4 

ls22s22p  (2P1/2) 

2 

2s22p3  (2D)(2P1/2) 

12 

2s22p  (2P3/2) 

4 

2s22p3  (2P3/2) 

4 

2s2p2  (4P) 

12 

2s2p4  (4P) 

12 

2s2p2  (2D)(2P1/2) 

12 

2s2p4  (2D) 

2s2p4  (2S1/2)(2P3/2) 

10 

6 

2s2P2  (2S)(2P3/2) 

2p3  C4S) 

6 

4 

2s2p4  (2P1/2) 

2 

2p3  (2D)(2P1/2) 

12 

2p5  (2P) 

6 

2P3  C2P3/2) 

4 

2sz2p^3£ 

270 

2s23£ 

18 

2s2p33£ 

720 

2s2p(l/2)3£ 

82 

2s22p24£ 

480 

2s2P(3/2)3 £ 

134 

2s2p34£ 

1280 

2s24£ 

32 

2s2p4£ 

384 

C-like 

Be-like 

ls22s22p2  (3Pq) 

1 

2s22p2  (3P1)(1D2) 

8 

1s22s2  (1Sq) 

1 

2s22p2  (3P2)(1Sq) 

6 

2s2p  (3P) 

9 

2s2p3  (5S2) 

5 

2s2p  (lP1) 

3 

2s2p3  (3D) 

15 

2p2  (3?q) 

1 

2s2p3(3PG  1)(1D)(3S) 

12 

2p2  (3P1)(1D2) 

8 

2s2p3  C3P2)(1P1) 

8 

2p2  (3P2)(1Sq) 

6 

2P4  (3P0>2) 

6 

2s3s/3p 

16 

2p4  (3P1)(1D2) 

8 

2s  3d 

20 

2p4  (LS0) 

1 

2p(l/2)3 £ 

44 

2s22p3£ 

108 

2p(3/2)3£ 

64 

2s2p23£ 

540 

2s4£ 

36 

2s22p4£ 

192 

2p4£ 

192 

2s2p24£ 

960 

14 


Table  I  —  Atomic  level  structure  used  in  the  ionization/radiation  model  for  neon, 
argon,  and  krypton  (only  the  krypton  model  has  Na-like  excited  levels). 


F-like 

— 

ls22s22p5  (2P) 

6 

2s2p6  ( 2S) 

2 

Na-like 

2s22p43s 

30 

2s22p43p 

90 

2s22p63s 

2 

2s22p43d 

150 

3p 

6 

2s2p53s 

24 

3d 

10 

2s2p33p 

72 

4s 

2 

2s2p53d 

120 

4p 

6 

2s22p44£ 

480 

4d 

10 

2s2p34£ 

384 

4f 

14 

2s22p45£ 

51 

50 

2s2p55£ 

1350 

61 

72 

2s22p46£ ) 

1944 

2s2pV  ) 

Ne-like 

— 

O-like 

2s22p6  (1S) 

1 

" 

2s22p53s  (3/2) 

8 

ls22s22p4  (3P2)(1S  ) 

6 

2s/2p:>3s  (1/2) 

4 

2s22p4  (3P1)(1D2) 

8 

2s^2p'>3p  (3/2) 

24 

2s22P4  (3Po) 

1 

2s/2p:>3p  (1/2) 

9  t; 

12 

2s2p5  (3P) 

9 

2s^2p^3d  (3/2) 

40 

2s2p5 

3 

2s^2p^3d  (1/2) 

A 

20 

2p6  (1SD) 

1 

2s2pb3s 

4 

2s22p33£ 

360 

2s2pb3p 

12 

2s2p33£ 

270 

2s2pb3d 

9  ? 

20 

2s22p34£ 

640 

2siJ2p;)4  z 

192 

2s2p44£ 

480 

2s^2p^5£ 

300 

2s22p56£ 

432 

13 
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Et  =  |  (1  +  Z)  Nt  kT 
T  2  I  e 

Because  the  average  Z  for  reaching  the  krypton  L-shell  (26)  is  greater  than 
that  for  argon  K-shell  (16),  more  energy  is  required  to  attain  the  same 
temperature.  Second,  some  of  the  kinetic  energy  is  converted  to  potential 
energy  via  ionization  of  the  plasma  ions  to  the  appropriate  electron 
shell.  To  reach  the  argon  K-shell  requires  only  about  5.9  keV/ion  of 
ionization  energy,  while  to  ionize  krypton  to  the  L-shell  requires  12.65 
keV/ion.  Third,  it  has  been  shown^1  that  to  effectively  "burn  through"  an 
electron  shell  is  not  a  simple  matter  of  supplying  enough  energy  for 
ionization.  Coilisional  excitations  and  radiative  decay  occur  on  a  time 
scale  several  orders  of  magnitude  faster  than  coilisional  ionization; 
hence,  radiation  loss  during  the  ionization  phase  must  be  taken  into 
account,  and  the  actual  energy  required  is  given  by 

E  =  EI  +  l  /  Pr  (Te(t))  dt  (6) 

where  Ej  is  ionization  energy  and  PRZ  is  the  radiation  loss  of  ion  z.  Just 
as  the  radiative  cooling  rate  for  the  krypton  L-shell  exceeds  that  of  the 
argon  K-shell,  (see  Fig.  1)  the  krypton  M- shell  cooling  rate  exceeds  that 
of  the  argon  L-shell.  Thus,  efficient  radiation  from  the  lower  shells  of 
higher-Z  materials  may  necessarily  require  more  powerful  generators  to 
attain  the  same  conditions  as  higher  shells  of  lower-Z  loads. 

Finally,  it  was  shown  that  "spectral  fill-in"  can  be  obtained  by 
mixtures  using  appropriate  concentrations  of  complementary  constituents, 
resulting  in  smoother  photon  distri butions  in  the  line  spectrum.  However, 
it  may  be  necessary  to  carefully  tailor  the  temperature  gradient  in  the 
dense  plasma  to  allow  all  the  relevant  ions  to  contribute  to  the  broadband 
emission  spectrum. 
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continuous  line  emission  from  0.9  to  4.0  keV.  In  fact,  at  plasma 
temperatures  below  600  eV,  the  argon  L-  and  krypton  M-shell  emission  would 
fill  in  the  spectrum  below  0.8  keV,  as  well;  it  would  take  very  careful 
tailoring  of  the  load  and  generator  characteristics ,  however,  to  maintain 
an  imploded  plasma  with  a  temperature  gradient  ranging  from  a  few  hundred 
electron  volts  up  to  2.5  keV! 

IV.  CONCLUSIONS 

An  investigation  of  the  radiative  yield  of  puffed-gas  mixtures  has 
been  initiated  in  which  the  temperature  and  concentration  effects  have  been 
studied.  The  main  result  of  the  simulations  is  that  the  yield  of  the 
mixtures  is  always  less  than  one  can  obtain  by  using  one  or  the  other  of 
the  constituents  in  pure  form.  Consequently,  the  only  hope  of  producing 
enhanced  emission  via  this  approach  would  be  if  the  plasma  mixture  attained 
a  higher  temperature,  higher  density,  or  increased  mass  participation  than 
that  obtainable  from  the  pure  constituents  (earlier  experiments  on  a  lower- 
power  machine  have  indicated  some  evidence  of  temperature  enhancement  using 
mixtures).  Although  that  study  is  the  next  step  in  our  investigation  of 
radiation  enhancement,  and  requires  a  combined  radiation/hydrodynamics 
model  to  characterize  the  plasma  dynamics,  some  general  statements  can  be 
made  regarding  the  radiative  performance  of  the  load. 

Essentially,  the  goal  of  maximizing  the  photon  yield  is  attained  by 
efficient  conversion  of  kinetic  energy  of  the  imploding  plasma  to 
radiation.  In  general,  the  hotter  and  denser  the  plasma  becomes  the 
greater  will  be  the  yield,  although  the  temperature  dependence  is  subject 
to  the  problem  of  burn-through  of  the  atomic  shell  structure.  This  is 
easily  overcome,  however,  by  going  to  higher  Z  materials.  In  fact,  it  has 
often  been  suggested  that  one  can  obtain  more  emission  in,  for  example,  the 
hv  >  1.5  keV  range  by  going  to  the  krypton  L-shell  as  opposed  to  the  argon 
K-shell .  From  Fig.  1  it  is  obvious  that  if  the  plasmas  can  be  made  to 
implode  to  the  same  temperatures,  densities,  and  sizes,  the  krypton  plasma 
radiates  at  a  much  higher  power  level  than  argon.  However,  to  obtain  the 
same  conditions  in  both  plasmas  requires  more  energy  for  krypton  since  a 
number  of  obstacles  must  be  overcome.  First,  some  of  the  plasma  kinetic 
energy  is  converted  to  thermal  energy. 


significantly  harder  in  the  mixture  than  in  the  pure  krypton  case. 
Presumably,  one  could  then  infer  that  much  higher  temperatures  were 


attained  in 

the 

mixture  than  in  pure 

krypton. 

Due 

to 

the  low 

concentration , 

the 

krypton  ion  density  was 

down  from 

3.5  X 

1019 

cm- ^  in 

pure  krypton  to  2.7  X  lO1^  in  the  mixture,  and,  hence,  so  was  the  total 
yield,  from  about  250  Joules  down  to  110  Joules  (most  of  the  photons  were 
of  hv  <  400  eV).  The  interesting  question  that  emerges  is  whether 

significant  temperature  increases  can  still  be  maintained  at  krypton 
concentrations  of  5*  or  greater.  As  seen  in  Fig.  6,  if  the  plasma 
temperature  is  below  1.0  keV,  an  increase  of  only  100  eV  can  greatly 
enhance  the  yield  of  photons  above  1.5  keV.  Thus,  smaller  generators  may 
be  able  to  enhance  the  production  of  krypton  L-shell  lines  by  using 
helium/ krypton  mixtures.  The  same  situation  would  apply  to  argon  K-lines 
(see  Fig.  5)  if  the  temperature  enhancement  was  evident  in  mixtures  of 
helium/ argon. 

c.  Spectral  Enhancement 


Another  advantage  of  using  gas-puff  loads  comprised  of  mixture  is 
the  so-called  "spectral  fill-in"  phenomenon.  By  clever  choice  of 
complementary  constituents,  the  natural  "gaps"  between  different  electron 
shell  emissions  that  appear  in  the  line-dominated  spectra  of  imploded  Z- 
pinch  plasmas  can  be  eliminated.  This  produces  a  smoother  distribution  of 
photon  energy  over  a  given  frequency  range  which  is  advantageous  for  a 
number  of  applications  of  the  radiation  source.  Note,  for  example,  the  way 
the  L-shell  emission  from  krypton  and  the  K-shell  emission  from  argon 
combine  in  Fig.  7  to  provide  almost  continuous  line  emission  from  1.6  keV 
to  4.0  keV.  The  spectrum  in  Fig.  7  was  generated  by  assuming  a  steady- 
state  temperature  profile  from  1.0  to  2.5  keV  exists  in  the  puffed-gas 
implosion.  This  is  a  necessary  requirement  in  order  to  obtain  radiation 
from  all  eight  ions  in  the  krypton  L-shell;  the  two  argon  K-shell  ions,  on 
the  other  hand,  radiate  over  a  wide  temperature  range  (see  Fig.  1). 
Similarly,  in  Fig.  8,  a  neon/krypton  mixture  can  be  used  to  obtain  line 
emission  from  0.9  keV  up  to  2.8  keV  with  the  appropriate  concentration  and 
temperature  conditions  (approximately  10:1  neon  to  krypton  at  about  1 
keV) .  Conceivably,  one  could  mix  neon,  argon,  and  krypton  to  obtain 


of  higher  temperature  or  higher  density,  or  more  mass  to  participate  in  the 
implosion  as  compared  to  the  pure  substances.  These  possibilities  will  be 
discussed  in  more  detail  in  the  next  section. 

Although  the  main  emphasis  here  has  been  on  photon  energies  above  a 
kilovolt,  for  the  sake  of  completeness,  the  total  yield  is  shown  in  Fig.  4 
for  neon/ argon  mixtures.  Again,  the  total  is  the  solid  lines  while  the 
broken  curves  represent  the  argon  emission.  Qualitatively,  Figs.  2  and  4 
and  very  similar  in  their  behavior.  Quantitatively,  we  see  that  the  lower 
temperatures  produce  more  argon  yield  than  the  higher  ones  up  to  2  keV, 
where  the  K-shell  emission  begins  to  peak  (see  Fig.  1).  This  is  due  simply 
to  the  fact  that,  at  this  density,  the  L-shell  of  argon  is  a  more  efficient 
radiator  than  the  K-shell  (see  Ref.  13  for  a  complete  discussion  of  L-  and 
K-  shell  competition  for  photon  production  as  a  function  of  density). 

The  last  two  mixtures  studied  were  helium/ argon  and  helium/ krypton . 
The  results  of  these  mixed  loads  are  shown  in  Figs.  5  and  6;  here,  the 
solid  lines  represent  the  total  radiation  from  the  higher-Z  constituent 
while  the  broken  lines  indicate  the  soft  x-ray  component  (no  helium 
emission  is  shown  since  that  constituent  is  fully  ionized).  Since  helium 
is  such  a  light  ion,  the  total  mass  of  the  mixture  is  almost  all  due  to  the 
higher-Z  material  at  even  rather  low  concentrations  (see  Table  II).  Thus, 
the  yield  curves  rise  quickly  in  Figs.  5  and  6  to  values  approaching  the 
maximum  at  concentrations  of  only  20% .  In  fact,  only  a  5%  concentration  of 
the  highei — Z  constituent  gives  a  yield  which  is  down  by  only  a  factor  of 
two  from  the  maximum.  Once  again,  the  same  message  is  communicated  by  the 
graphs:  no  radiation  enhancement  is  possible  over  that  of  the  pure  gas 
load  by  using  a  helium  mixture  unless  the  lower  concentrations  can  be 
imploded  to  higher  density  or  higher  temperature. 

In  an  experiment  on  a  smaller  generator  (5  kj  at  30  kV),  Bailey1"^ 
reported  on  yield  measurements  for  pure  helium,  helium/ krypton  mixtures, 
and  pure  krypton  where  the  ion  number  concentration  of  krypton  was  kept  at 
1.5$  for  the  mixture.  The  most  significant  observation  was  a  large 
increase  in  the  total  emission  over  that  of  pure  helium.  Measured  electron 
densities  were  about  1020  cm“^,  five  times  that  of  pure  helium,  and  the 
minimum  radius  about  a  factor  of  two  smaller.  Also,  it  was  determined  that 
the  average  ionization  stage  was  higher  and  the  radiation  spectrum 


since  the 
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yields  shown  in  Fig.  1  scale  with  density  approximately  as  Mj 
relevant  ion  ground  states  are  close  to  coronal  equilibrium  at  the 
densities  commonly  encountered  in  gas-puff  Z-pinches.  Note  that  while  the 
argon  X-sheli  and  krypton  L-shell  line  emission  are  both  near  their  maxima 
in  the  temperature  range  of  2-6  keV,  the  neon  and  argon  will  not  emit 
strongly  together  at  the  same  temperature  in  the  photon  energy  region  of 
interest,  intimating  that  argon/ krypton  mixtures  may  be  better  radiators. 
On  the  other  hand,  based  on  past  performance  of  the  large  pulse  power 
generators,  it  is  not  clear  that  plasma  loads  can  be  heated  in  bulk  to 
temperatures  much  above  2  keV.  Thus,  actual  radiative  performance  of  the 
two  mixtures  may  not  differ  significantly  in  actual  machine  conditions. 

Considering  now  the  radiative  output  of  the  mixtures,  the  results  are 
presented  as  the  yield  (in  the  relevant  photon  range)  as  a  f 'unction  of  ion 
number  concentration  of  the  higher-Z  constituent  for  several  temperatures 
spanning  the  operational  range  of  the  machines.  The  results  of  the 
neon/ argon  study  are  shown  in  Fig.  2  for  hv  >  0.9  keV;  the  solid  lines  are 
the  neon  plus  argon  output  while  the  output  of  argon  alone  i3  depicted  by 
the  broken  lines  (thus,  both  curves  at  the  same  temperature  merge  at  1 00% 
argon  concentration) .  In  Fig.  2,  it  is  shown  that  the  neon  radiation  falls 
off  above  500  eV;  this  explains  the  behavior  of  the  yield  with  temperature 
at  0%  concentration  of  argon.  Conversely,  argon  yield  rises  with 
temperature  above  500  eV,  explaining  the  cross-over  of  the  solid  curves. 
As  expected  the  argon  emission  rises  rapidly  with  concentration  and 
asymptotes  with  the  neon  +  argon  curve.  The  main  point  to  be  emphasized 
from  Fig.  2,  however,  is  that  no  possible  mixture  of  the  two  constituents 
results  in  a  soft  x-ray  yield  which  exceeds  the  larger  of  the  two  pure 
yields.  Stating  it  another  way,  under  the  assumption  that  all  loads  of  a 
fixed  mass  of  two  constituents  will  implode  to  the  same  temperatures,  there 
appears  to  be  no  way  to  enhance  the  yield  of  a  load  by  using  mixtures  above 
that  wnich  could  be  obtained  from  one  or  the  other  of  the  pure 

substances.  This  conclusion  is  verified  in  Fig.  3  for  argon/krypton 

mixtures  (with  hv  >  1.5  keV) ,  where  very  similar  behavior  of  the  yield  with 
concentration  is  seen.  This  is  a  rather  unexpected  result  from  the  point 
of  view  of  radiation  enhancement  of  super  kilovolt  photons  using  mixed 
loads.  Of  course,  the  gas  mixture  may  still  be  shown  to  affect  the 
implosion  dynamics  in  a  positive  way,  such  as  allowing  either  an  attainment 
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densities  of  2-5  x  10  7  are  not  uncommon  for  high-yield  shots.  Having 
fixed  the  total  mass  density,  the  actual  ion  densities  used  in  the 
calculation  are  given  in  Table  II  for  the  concentrations  analyzed. 

At  this  point,  it  may  be  worthwhile  to  point  out  that  this  strategy  of 
parameter  variation  may  differ  from  actual  experimental  procedures.  In  the 
laboratory  it  is  common  to  obtain  cylinders  of  gas  mixtures  at  fixed 
concentrations ,  and,  by  varying  the  diaphragm  diameter  on  the  fast-opening 
valve  of  the  gas-puff  nozzle,  study  the  radiation  yield  as  a  function  of 
initial  mass  density.  The  reason  for  not  adopting  this  technique  is  due  to 
our  stationary  equilibrium  model.  Varying  the  initial  mass  of  the  load 
will  likely  alter  the  implosion  dynamics  and,  hence,  the  final  temperature 
and  density.  Thus,  a  meaningful  calculation  of  the  load  performance  would 
require  complete  modeling  of  the  flow  effects  which  is  beyond  the  intent  of 
this  paper.  However,  we  plan  to  eventually  couple  the  radiation/ 
ionization/mixture  model  to  a  1-D  MHD  code  for  the  purpose  of  assessing  the 
cooperative  radiati ve-hydrodynamic  effects  on  radiation  enhancement. 

b.  Radiative  Yields 

As  was  discussed  in  the  previous  section,  the  radiation  output  of 
the  mixtures  was  calculated  as  a  function  of  the  electron  temperature  and 
constituent  concentration  for  fixed  mass  density.  Two  radiative  yield 
measurements  are  of  interest:  the  total  radiative  output,  and  that  above  a 
spectral  cutoff  energy,  900  eV  for  neon,  1.6  keV  for  krypton,  and  3.0  keV 
for  argon.  Thus,  the  yield  above  cutoff  will  be  the  spectral  output  at 
about  a  kilovolt  and  greater,  since  no  strong  krypton  lines  occur  between 
900  and  1600  eV  and  no  strong  argon  lines  occur  between  900  and  3000  eV. 
Studies  of  radiation  enhancement  via  mixtures  at  photon  energies  below  1.0 
keV  require  different  experimental  configurations  and  are  not  considered  in 
optimum  conditions  here.  In  order  to  acquaint  the  reader  with  the 
temperature  variation  of  radiative  yields  for  pure  gases,  the  optically 
thin  power  densities  for  (a)  neon  K-lines  (hv  >  0.9  keV),  (b)  argon  K- 
lines  (hv  >  3.0  keV),  and  (c)  krypton  L-lines  (hv  >  1.6  keV)  are  shown  in 
Fig.  1  for  a  fixed  density  of  101®  ions/cm^.  This  figure  will  be  very 
useful  in  explaining  the  variation  in  yields  from  the  mixture  gases  as  the 
electron  temperature  of  the  mixture  is  increased.  The  line  radiation 


7 


%  CONCENTRATION  ARGON 

Fig.  5.  Radiation  yield  versus  argon  concentration  for  helium/argon 

mixtures;  solid  curves  are  total  radiation,  broken  curves  are  for 
hv  >  3-0  keV. 


HELIUM/ KRYPTON 


Radiation  yield  versus  krypton  concentration  for  helium/krypton 
mixtues;  solid  curves  are  total  radiation,  broken  curves  are  for 


KRYPTON/ARGON 


PHOTON  ENERGY  (keV) 

Fig.  7.  Radiation  spectrum  resulting  from  an  imploded  argon/ krypton 
mi xture . 


NEON/KRYPTON 


Fig.  8.  Radiation  spectrum  resulting  From  an  imploded  neon/krypton 
mi xture . 
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B.  Modeling  of  K-Shell  Yields 


I.  INTRODUCTION 

It  nas  long  been  a  DNA  goal  to  maximize  the  K-shell  x-nay  yield 

obtained  in  plasma  discharges  created  with  pulsed-power  generators  such  as 

BLACKJACK  5  and  PITHON.  During  FY85  important  strides  were  taken  in  both 

comprenending  previous  results  and  in  new  concept  development. 

A  key  achievement  in  assessing  the  contributory  factors  to  previously 
obtained  K-shell  yields  and  projecting  the  yields  expected  from  future, 
higher-power  machines,  is  described  in  the  attached  NRL  Memorandum  Report 
No.  5406,  "K-Shell  Yield  Scaling  Law  for  Conventional  PRS  Loads."  This 
report  is  self-explanatory  and  requires  no  elaboration  here. 

We  have  also  performed  a  detailed,  quantitative  assessment  of  the 

possibility  of  yield  enhancement  through  use  of  optically  thick  mixtures  of 
elements  with  similar  atomic  number,  and  find  the  concept  to  be  very 

promising.  In  a  similar  study  of  optically  thin  mixtures,  discussed 

elsewhere  in  this  report,  it  was  found  that  no  yield  increase  was  possible, 
but  the  spectral  distribution  could  be  tailored  to  the  needs  of  the 

simulation  by  employing  elements  of  widely  divergent  atomic  number. 

The  present  study  is  concerned  with  a  considerably  different  mixture 
scenario:  the  plasma  is  optically  thick  and  the  atomic  numbers  very 

close.  In  this  case,  opacity  effects  do  influence  the  total  yield  at  the 
expense  of  spectral  filling  since  the  spectral  lines  from  the  various 
shells  of  the  two  elements  lie  fairly  close  in  energy.  As  seen  below,  the 
effects  produced  by  this  type  of  mixture  differ  considerably  from  optically 
tnin  mixtures  of  significantly  different  atomic  number. 

We  consider  sodium  and  neon  ions,  mixed  in  varying  ratio,  and  heated 
to  temperatures  at  which  the  helium  and  hydrogen-li ke  ionization  stages  are 
by  far  the  dominant  species.  The  choice  of  sodium  and  neon  is  of 
significant  interest  not  only  because  of  the  close  atomic  numbers,  but  also 

p  i 

because  the  Na  X  Is  -ls2p  P  line  at  11.0027  A  virtually  coincides  with 
Ne  IX  1 s 1 3 4p ^  p  at  1  1.0003  A.  Therefore,  the  effect  of  direct  overlap  of 
high-opacity  spectral  regions  may  be  ascertained. 
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The  temperatures  and  densities  attained  in  a  pinch  with  a  given 
current-voltage  pulse  are  largely  dependent  on  the  plasma  mass  which  is 
accelerated.  Therefore,  a  reasonable  comparison  of  the  radiative 
emissivity  of  a  mixture  with  that  of  a  one-element  plasma  may  be  made  by 
holding  the  mass  density  constant.  Since  sodium  and  neon  have  nearly  the 
same  mass,  in  the  cases  presented  below  we  have  kept  the  total  ion  density 
constant  while  varying  the  sodium/neon  ratio,  to  remain  reasonably  faithful 
to  the  prescription  for  a  meaningful  comparison. 

For  a  plasma  of  a  given  3ize,  total  ion  density,  temperature  and 
sodium/neon  ratio,  we  solve  for  the  level  populations  and  emitted 
radiation.  The  plasmas,  which  persist  for  tens  of  nanoseconds,  are  assumed 
to  be  in  collisional-radiation  equilibrium  (CRE),  that  is,  the  level 
populations  reach  values  which  permit  them  to  be  in  a  steady-state 
consistent  with  the  other  level  populations,  atomic  rates,  and  the  ambient 
radiation  field.  To  perform  this  calculation  the  critical  elements 
required  are  the  atomic  level  model  and  rates,  and  a  radiation  transport 
algorithm  which  accounts  for  two-element  opacity  overlap. 

We  consider  only  sodium-neon  plasmas  of  such  temperatures  and 

densities  that  the  K -shell  dominates  the  ionic  species.  Hence,  only  ground 

states  are  carried  in  the  model  for  all  stages  below  lithium-like.  For  the 

lithium-like  stages,  the  1  s22s  2S  (ground),  1s22p  2P,  1s23s  2S,  1s23p  2P, 
P  P  2  P 

Is  3d  D,  Is  and  Is  51  levels  are  included.  The  helium-like  stage 

T  1  3 

carries,  in  addition  to  the  ground  state,  the  1s23  JS,  1s2s  S,  1s2p  JP , 
1  s2p  1  P ,  1s3H  ,  1s3i  1L,  1 3 4 2, ,  and  135£  levels.  The  hydrogen-like  levels 
include  n  =  1  through  n  =*  5 . 

The  atomic  processes  populating  and  depopulating  the  levels  are: 
spontaneous  radiative  decay,  electron  collisional  excitation  and 
deexcitation,  collisional  ionization,  3-body,  radiative,  and  dielectronic 
recombination.  The  continuum  is  optically  thin  in  these  plasmas;  hence 
photoionization  is  neglected.  Line  photoexcitation  is  calculated  by  a 
self-consistent  radiative  transfer  model. 
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a.  K- Shell  Output  of  Optically  Thick  Mixtures 


Calculations  have  been  performed  for  a  variety  of  assumed  plasma 

conditions  to  determine  the  total  K-shell  resonance  line  power  as  a 

function  of  the  sodium/neon  ratio.  In  each  case  a  cylindrical  plasma  of 
length  4  cm  of  fixed  total  (sodium  plus  neon)  ion  density  is  assumed.  The 
temperature,  diameter,  and  sodium/neon  ratio  were  varied,  and  the  results 
are  summarized  in  Fig.  1  and  Table  I.  In  Fig.  1  the  K-shell  power  output 
in  J/nsec  is  plotted  as  a  function  of  the  fractional  abundance  of  sodium 
ions  for  a  variety  of  plasma  conditions.  Note  that  generally  the  power 
output  (equivalent  to  K-shell  emissivity)  peaks  at  a  sodium  ion 
concentration  of  50-75?,  and  that  the  emissivity  in  some  cases 

substantially  exceeds  that  obtained  either  for  pure  neon  or  for  pure 

sodium.  The  effect  is  more  pronounced  as  the  line  opacity  of  the  plasma 
(roughly  proportional  to  the  product  of  ion  density  and  diameter) 

increases.  Therefore,  the  K-shell  emissivity  of  an  opaque  plasma  mixture 

may  exceed  that  for  either  element  alone  for  similar  atomic  number 

elements.  The  precise  nature  of  this  effect  is  illuminated  in  Table  I, 

where  the  power  outputs  for  sodium  and  neon  Ly  a  and  He  a  (Is  -1s2p  P) 
lines  are  given  for  various  plasma  conditions.  These  two  lines  usually 
account  for  more  than  65?  of  the  K-shell  power.  Were  the  composition  of  an 
optically  thin  plasma  to  be  changed  from  pure  neon  to  50?  neon,  for 

instance,  the  line  power  output  would  decrease  by  50?  -  ignoring  small 
changes  in  the  electron  density  due  to  sodium.  However,  inspection  of  the 
table  reveals  that  the  neon  lines'  output  never  drops  by  as  much  as  50?  and 
in  fact  usually  drops  by  only  about  20?.  The  one  exception  is  the  least 
opaque  plasma  of  ion  density  lO1^  cm"^,  where  the  Ne  IX  He  a  power  drops  by 
nearly  40?.  In  this  least  opaque  plasma  the  Ly  a  optical  depth  is  -50. 
The  same  effect  exists  for  sodium.  The  change  of  the  photon  output  is  due 
to  opacity  effects.  This  results  in  the  greatest  emissivity  occuring  for 
an  intermediate  mixing  ratio  of  sodium-to-neon,  where  there  are  twice  as 
many  lines  radiating  as  in  the  pure  element  case  and  the  optical  depths  of 
the  radiating  lines  are  smaller  than  for  the  pure-element  case. 
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b. 


Radiation  Physics  Controlling  Power  Output 


In  analyzing  the  transfer  of  line  photons  in  laboratory  plasmas,  the 
following  conceptual  picture  is  often  useful.  A  line  photon  is  "created" 
when  the  upper  level  cf  the  transition  is  collisionally  populated  and  then 
decays  radiatively.  Following  its  creation  in  an  optically  thick  plasma 
the  photon  will  generally  be  reabsorbed  and  reemitted  a  considerable  number 
of  times  before  being  quenched  or  ultimately  escaping  the  plasma.  During 
each  reacsorotion  there  is  a  finite  probability  P^  that  the  level  will  be 
collisionally  depopulated  and  t.oe  line  photon  destroyed.  This  probability 
?q  is  given  by  D/(D+A)  where  DCsec*1)  is  the  sum  of  the  coliisional  rates 
out  of  the  upper  transition  level  times  the  electron  density  and  A  is  the 
spontaneous  decay  rate  for  the  line.  Following  its  initial  coliisional 
creation  t.oe  photon  will  jltimately  suffer  one  of  two  fates:  destruction 
or  eventual  escape  from  the  plasma.  *e  denote  the  probaDility  of  eventual 
escape  oy  the  probability  of  escape  on  any  single  flight  is  given  by 
PQ.  The  quantity  Pu  may  be  obtained  by  considering  the  condition  for  a 
steady-state  population  N  of  t.oe  upper  level.  Letting  C  stand  for  the 
coliisional  creation  rate  this  condition  is 

C  *  N  ( D+AP  )  (cm3/sec)  (1  ) 
e 


where  in  Eq.  (1),  radiative  reabsorption  of  the  line  photon  is  taken  into 
account  by  diluting  the  decay  rate  by  the  single-flight  photon  escape 
probability.  Even  in  a  uniform  plasma,  Pg  varies  with  position  -  Eq.  (1) 
is  intended  as  spatially  averaged  approximate  description  of  the  photon 
transfer.  After  a  few  simple  manipulations,  Eq.  (1)  yields,  taking  into 
account  that  Pn  =  D/(D+A), 
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However,  the  left-hand-side  of  Eq.  (2)  is  the  total  photon  emission  rate 
divided  by  the  total  upper  level  excitation  rate;  that  is,  precisely  the 
ultimate  escape  probability.  Therefore, 
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in  this  approximation. 

In  Figs.  (2)  and  (3)  we  have  plotted  the  quantities  ar.d 

Pa(  1  -Pq ) /'[Pq  +  ( 1  -Pq)Pa]  vs.  sodium  ion  density  and  diameter,  respectively. 

These  parameters  are  displayed  for  the  Na  X  is^-isZp?  line  for  a 

temperature  of  400  eV  and,  for  Fig.  2,  a  fixed  diameter  of  1.4  mm.  Fig.  3 

TO 

reflects  calculations  f.r  a  fixed  sodium  ion  density  cf  lO*-'"'  cm  .  ?a  and 
?r>  are  obtained  from  Voigt  profile  formulae  and  the  atomic  rates,  respec- 
tively.  The  accuracy  of  the  spatial-average  approximation  is  indicated  in 
both  Figs.  2  and  3  by  also  plotting  P  as  numerically  calculated  from  the 
actual  computed  line  photon  output  power.  As  is  obvious,  excellent 

correspondence  is  obtained  with  the  probabilistic,  spatial-average 
approximation  of  Eq.  (3)*  This  is  especially  significant  in  light  of  the 
fact  that  even  though  the  plasma  temperature  and  density  are  assumed 
spatially  uniform,  the  presence  of  a  boundary  guarantees  a  signficant 
degree  of  nonuniformity  in  both  the  photon  field  and  the  spatial  profile  of 
the  level  populations. 

The  important  feature  of  Fig.  3  relevant  to  the  mixture  emissivity  is 

the  behavior  of  Pu  with  respect  to  opacity.  Optical  depth  varies 

approximately  linearly  with  plasma  diameter.  The  electron  density  in  a 

sodium-neon  plasma  is  very  nearly  constant  with  sodium/neon  mixing  ratio 

provided  the  total  ion  density  is  constant.  Therefore,  halving  the 

diameter  of  a  pure  sodium  plasma  would  lead  to  an  increase  in  Pu  very 

similar  to  that  which  is  obtained  by  comparing  a  pure  sodium  plasma  to  one 

of  the  same  ion  density  which  is  half  sodium,  half  neon.  The  profile  of 

-1  /2 

Fig.  3  indicates  that  Pp  ~  t  for  optically  thick  K-shell  sodium 

plasmas  of  ion  density  10  cm  J  in  the  1-10  mm  range.  Therefore,  halving 
the  optical  depth  by  reducing  the  sodium  fraction  from  1.0  to  0.5  results 
in  Pu  increasing  by  401 .  The  other  sodium  and  neon  lines  experience  a 
similar  -  through  not  identical  -  increase  in  Pu,  which  accounts  for  the 
curves  of  Fig.  1 . 
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replications  for  DMA  Objectives 


The  above  results  show  convincingly  that  the  K-sheli  emissivity  of  a 
plasma  may  be  increased  by  employing  mixtures  of  elements  of  similar  atomic 
number  -  where  the  total  mass  employed  in  a  machine  such  as  BLACKJACK  5 
•would  remain  the  same  as  that  used  on  a  conventional  single-element  shot. 
The  mechanism  that  increases  the  emissivity  is  that  the  optical  depths  of 
each  element  are  less ,  allowing  more  photons  out  -  and  twi ce  as  many  lines 
are  present  in  the  K-shell  region.  The  choice  of  sodium  and  neon  does  not 
mean  that  these  promising  results  would  be  confined  to  these  elements. 
Indeed,  the  physics  would  be  very  similar  for  mixtures  such  as  aluminum- 
silicon,  magnesium- aluminum,  argon-titanium,  or  calcium-titanium.  We 
recommend  that  a  series  of  experiments  using  one  or  more  of  the  above 
mixtures  be  carried  out  to  test  these  concepts.  Conceivably,  the  use  of 
three  ;r  more  elements  could,  by  further  enhancing  the  photon  escape  ana 
multiplying  the  number  of  K-shell  lines,  result  in  even  greater  increases 
in  power  output  than  calculated  here  for  the  sodium/neon  mixture. 
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TABLE  I:  Power  (Outputs  (J/nsec)  for  Lv  a  and  He  a 


d  =  5.6  mm,  Te  =  400eV,  Mj  =  5xl0^9  cm  3 
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XI 
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a 
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— 

d  =  1.4 

mm,  Tg  =  400 eV 

,  Nj  =  1020 

cm-3 

Pure  Ma 

75Na, 25Ne 

50Na, 50Ne 

75Ne,  25Na 

Pure  1 

Ne 

IX 

He 

a 

— 

14. 
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K  SHELL 
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5  x  102 


Fig.  1.  The  total  K-shell  power  output  in  J/nsec  from  a  sodium-neon  plasma 
is  plotted  against  fractional  abundance  of  sodium  for  a  range  of 
indicated  plasma  conditions.  In  each  case  the  plasma  is  assumed  to  be 
a  4  cm-long  cylinder  with  fixed  total  ion  density  (sodium  plus  neon 
ions) . 


A  typical 

problem  mesh  consists  of  32  z 

ones  in  the 

R  direction 

and  1  c 

zones  along  2. 

Tlais 

relatively  coarse  grid 

is  dictated 

at  present 

by  both 

memory  and  run 

time 

considerations  of  the 

atomic  physics  model. 

After 

initialization,  the  Z-D  PRISM  node  was  run  witn  the  atomic  physics  pac' age 
timed  on.  Although,  at  present,  there  is  no  radiation  transport  in  the  2- 
0  runs.  As  a  comparison,  several  runs  were  also  made  without  t.ne  atomic 
physics  model. 

Figure  3  shows  density  contours  out  to  the  95?  mass  line  (i.e.,  95?  of 

the  mass  lies  inside  of  the  last  contour)  at  1.12  usee.  The  2-D  run  was 

initialized  at  i .05  jsec.  In  this  figure  and  in  the  figures  which  follow, 

the  wavelength,  \,  is  1  cm.  The  contours  are  separated  logari thmi tally  by 

0.25  or  cy  one  quarter  of  a  decade.  As  discussed  above,  the  shell  is 

roughly  1  cm  in  thickness  and  can  be  defined  by  one  full  decade  (i.e.,  four 

contour  lines)  centered  near  3.5  cm.  The  backside  has  begun  to  show  signs 

of  an  instability  at  1.13  ysec  as  evidenced  by  the  contours  near  R  =  3  cm 

(2=0)  in  Fig.  4.  3y  1.24  usee,  the  disturbance  has  penetrated  to  the 

fourth  countour  as  shown  in  Fig.  5.  Figure  6  shows  density  contours 

computed  from  the  mesh  at  the  same  time  but  the  contours  only  go  out  to  the 

30?  mass  line.  The  same  quarter  decade  contour  separation  in  the  density 

contours  has  been  used  so  that  comparisons  can  be  made.  It  can  be  seen 

that  at  this  time  the  region  near  the  30?  mass  line  has  only  been  slightly 

perturbed.  In  Fig.  7,  the  contours  are  shown  just  prior  to  assembly  on 

axis.  Contours  are  again  snown  out  to  the  95?  mass  line  and  indicate  that 

the  perturbation  has  grown  significantly  at  the  rear  of  the  shell. 

However,  the  interior  portion  of  the  shell,  which  at  this  time  lies  within 

1  cm,  still  remains  unperturbed  and  intact.  Figure  8  shows  the  contours 

out  to  the  30?  mass  line  at  this  time.  At  this  time,  the  disturbance  has 

penetrated  only  slightly  more  into  the  shell  than  at  1.24  psec.  We  would 

expect  this  from  looking  at  the  1-0  plots  in  Fig.  2  where  it  is  seen  that 

even  near  the  time  of  implosion,  the  regions  where  rB,  and  p  cross  is  verv 

0 

narrow.  This  is  characteristic  of  a  "snowplow"  implosion  where  all  the 
mass  of  the  snell  3nd  tne  background  gas  interior  to  the  shell  is  swept  up 
by  the  field.  Since  this  effect  creates  only  a  narrow  region  with  large 
gradients  near  the  backside  of  the  shell,  this  has  been  referred  to  as 
"snowplow  stabilization".  In  comparison,  Fig.  9  shows  the  shell  out  to  the 
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b.  Two-Dimensional  Implosion  Modeling 

An  important  issue  for  nigh  velocity  Z-pinch  PRS  loads  is  the 
stability  of  the  plasma  shell  as  it  is  accelerated  under  the  action  of  the 
axiai  magnetic  field.  In  the  case  of  an  accelerating  shell,  the  backside 
is  unstable  in  the  region  where  the  density  gradient  is  of  opposite  sign  to 
the  acceleration.  The  issues  are  then  twofold:  (1)  how  much  of  the  shell 
mass  is  in  this  unstable  configuration;  and  (2)  what  is  the  magnitude  of 
the  growth  rate  of  the  instability  ana  how  does  this  compare  to  the 
implosion  time.  A  summary  of  our  results  obtained  from  modeling  a 
representati ve  load  is  given  below.  A  more  detailed  description  will 
follow  under  separate  cover. 

The  load  we  have  chosen  to  model  is  a  120  pgm/cm  argon  gas  shell 
initially  at  20  cm  and  with  a  width  of  roughly  3  cm.  Ten  percent  of  the 
mass  lies  within  15  cm  and  represents  the  background  gas.  The  current 
pulse  rises  to  3*6  MA  at  1  psec  and  then  drops  to  a  small  value  by  1.3 
usee . 

We  modeled  this  problem  in  two  steps.  First,  a  1-D  PRISM  run  was  made 

with  all  the  relevant  physics,  including  radiation  transport,  turned  on. 

Second,  the  output  from  this  run  was  then  used  to  initialize  several  2-D 

runs.  Figure  1  shows  a  plot  of  several  mass  radii  versus  time.  The  mass 

percentile  values  indicate  how  much  mass  lies  inside  that  mesh  line.  From 

this  plot,  we  can  see  that  the  shell  does  not  become  fully  compressed  and 

accelerated  'until  about  0.9  usee .  This  is  slightly  before  peak  current  . 

Figure  2  shows  density  and  r B.  profiles  for  three  separate  times.  At  1.05 

9 

usee  tne  width  of  half-maximum  of  the  shell  is  approximately  1  cm .  A 
density  "foot"  is  seen  to  extend  another  1  cm  in  front  of  the  shell.  The 
backside  region  where  the  gradients  are  in  opposite  directions  lies  within 
a  region  less  than  1  cm. 

The  2-D  runs  were  initialized  by  expanding  the  1-D  profiles  at  a 
selected  time  in  the  Z-direction  and  then  adding  a  perturbation  which  is 
centered  around  the  density  peak  of  the  imploding  shell.  Symmetry  allows 
tne  Z  domain  to  extend  from  0  to  ,\/2  where  A  is  the  perturbation 
wr/olength.  The  quantity  A/\,  where  A  is  the  amplitude  of  the 
pert ur bat  ion ,  was  fixed  at  0.1  for  all  wavelengths  considered. 
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Fig.  5  —  GAMBLE  II  implosion  time  as  a  function  of  mass . 
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TIME  Crt5ec) 

Fig.  3  —  Time  history  of  several  mass  lines  for  a  20  ngm/cm  GAMBLE  II  implosion. 
Also  shown  is  the  time  history  of  the  total  current,  I. 
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Fig.  2  —  Density  and  magnetic  field  profiles  at  two  separate  times 
for  a  20  /ugm/cm  GAMBLE  implosion. 


43 


I  4.0 1  i  I  ^  F?HO 

1  i 


0.2  0.4  0.0  0.8  1.0  1.2 


R 

Fig.  1  —  Density  and  temperature  profiles  at  two  separate  times  for  a  20  ngm! 
cm  GAMBLE  II  implosion.  The  cross-matching  indicates  the  region  where 
radiation  is  being  absorbed. 
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2  kj  over  the  4  cm  length.  Overall,  the  yield  values  displayed  in  this 
figure  overlap  those  given  by  the  experiment.  The  qualitative  picture  i3 
of  a  peaking  of  the  yield  with  mass,  followed  by  a  rollover  of  the  yield 
curve.  At  low  mass,  the  yield  is  a  function  of  the  amount  of  material 
which  is  available  to  radiate.  At  .oigher  masses  (i.e.  greater  than  30 
ugm/cm  it  becomes  increasingly  more  difficult  to  accelerate  the  shell  so 
that  the  mass  is  driven  to  lower  peak  temperatures  during  the  implosion 
process . 


The  dispersion  near  the  peak  yield  is  due  to  the  fact  that  at  an 

optimum  temperature,  the  yield  is  sensitive  to  both  the  density  at  that 

temperature  as  well  as  the  amount  of  mass.  Since  the  power  goes  as  the 

density  squared,  the  integration  over  an  effective  volume  and  time  interval 
gives  a  yield  which  is  proportional  to  the  average  density.  Thus,  changes 
in  the  average  density  of  roughly  20 %  can  account  for  the  difference 
between  the  mean  yield  and  the  maximum  yield  if  the  temperature  is  assured 
to  be  constant  during  the  implosion.  It  should  be  remembered  that  the 
density  changes  by  roughly  two  orders  of  magnitude  during  the  implosion. 
Also,  we  have  ignored  important  variations  due  to  temperature,  opacity,  and 
other  dynamical  effects  in  the  above  discussion.  Therefore,  the  best  that 
one  can  hope  for  out  of  a  study  such  as  this  is  a  qualitative  picture  which 
quantifies  values  within  a  sensitivity  range  as  shown  in  Fig.  4. 

Figure  5  shows  the  implosion  time  as  a  function  of  mass.  Again,  a 

large  number  of  runs  were  used  in  making  this  plot  and  a  spread  of 
approximately  10  nsec  is  seen  for  a  given  mass  value.  The  implosion  times 
for  masses  in  the  20  ugm  to  35  ugni  range  are  consistent  with  the  values 
obtained  in  the  experiment. 


The  time  history  of  the  implosion  for  several  mass  radii  is  shown  in 
Fig.  3.  The  current  waveform  is  also  snown.  A  snowplow  is  seen  to  develop 
between  40  r.sea  and  60  nsec.  It  is  important  to  note  that  the  plasma  is 
still  at  902  of  its  initial  radius  at  the  time  of  peak  current.  The 
implosion  itself  occurs  at  110  nsec  when  the  total  current  is  substantially 
reduced.  Disassembly  occurs  rapidly  after  this  time.  The  shell  is  within 
302  of  its  initial  radius  for  30  nsec  and  this  roughly  defines  the  time 
span  of  the  radiation  pulse. 

Figure  9  shows  the  power  output  as  a  function  of  time.  Four  curves 
are  shown  and  represent  total  yield  and  K-line  yield  for  the  20  ygm/cm  case 
with  the  radiation  transport  turned  on  and  with  the  radiation  transport 
turned  off.  There  are  several  salient  features  in  this  figure.  First,  one 
notes  that  in  the  case  without  transport,  there  is  significant  power  output 
between  9q  nsec  and  90  nsec.  This  feature  dissappears  when  the  radiation 
transport  is  turned  on.  Most  of  this  radiation  comes  from  the  L-shell  but 
is  severally  inhibited  by  opacity  effects  in  the  case  with  transport.  This 
is  evidenced  by  the  previously  discussed  cross-hatched  region  that  appears 
in  Fig.  1 . 

Another  feature  is  the  difference  in  the  implosion  time  and  peak  power 
output  between  the  two  cases.  Again,  these  effects  can  be  explained  by  the 
fact  that  the  opacity  tends  to  keep  the  temperature  higher.  This  tends  to 
slow  down  the  implosion  because  of  an  increased  ''back"  pressure  pushing 
against  the  inward  moving  gas  shell.  The  larger  K-shell  output  occurs 

because  the  fraction  of  the  shell  which  contributes  to  the  radiative  yield 
during  the  implosion  phase  reaches  a  higher  temperature  with  opacity 

effects  included. 

Figures  9  and  5  show  the  composite  results  of  roughly  25  computer 
runs.  Variations  were  made  in  mass  (increasing  along  the  horizontal  axis 
in  both  figures)  as  well  as  turning  the  radiation  transport  on  and  off, 
adjusting  a  multiplier  on  the  magnetic  field  diffusion  coefficient,  and 

varying  the  initial  temperature  of  the  plasma. 

The  K-iine  yield  displayed  in  Fig.  9  shows  a  peak  near  22  ugm/cm  with 
a  mean  value  of  950  j/cm.  Over  a  9  cm  experimental  length,  this 

corresponds  to  a  total  yield  of  1.8  kJ.  By  varying  the  parameters  in  the 
proper  fasnion,  it  can  be  seen  that  it  is  possible  to  obtain  greater  than 
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C .  Dynamics  of  Imploding  Z-Pinch 


I.  GAMBLE  II  MODELING 

During  the  past  year,  we  have  modeled  the  Gamble  II  experiment  with  a 

one-dimensional  version  of  our  radiation  hydrodynamics  simulation  code, 

PRISM.  This  version  of  the  code  employs:  (1)  radiation  transport  where 

opacity  effects  are  important,  (2)  equation  of  state  utilizing  non-LTE 

atomic  physics  and  (3)  an  input  current  wave  form  coupled  with  a  diffusive 

equation  for  the  magnetic  field,  3  .  We  have  also  used  the  two- 

0 

dimensional  version  of  PRISM  to  model  the  later  pnase  of  tne  implosion  near 
axis.  At  present,  this  version  of  the  code  does  not  transport  radiation 
and  is  therefore  only  valid  in  the  thin  limit. 

a.  1-D  Modeling 

The  NRL  Gamble  II  neon  gas  puff  experiment  has  a  current  risetime  of 
approximately  65  nsec  with  a  peak  current  of  1.25  MA.  The  initial  central 
radius  of  the  gas  shell  formed  by  the  supersonic  nozzle  is  assumed  to  be 
1.25  cm  with  a  3  mm  width.  Thus,  a  mass  of  20  ygm/cm  corresponds  to  an 
initial  density  of  7.31  ygm/cm^  in  the  shell.  Here  we  have  assumed  that 
the  density  inside  the  shell  is  ten  percent  that  of  the  shell. 

Figure  1  shows  the  evolution  of  a  Z-pinch  implosion  utilizing  Gamble 

II  parameters.  Density  and  temperature  profiles  are  displayed.  Two  times 

are  shown  corresponding  to  a  time  near  peak  current  and  a  time  near 

implosion.  The  cross-hatching  under  the  temperature  profile  indicates  that 

this  region  is  a  net  absorber  of  radiation.  The  effect  of  this  feature  is 

to  help  maintain  the  temperature  of  the  core.  Figure  2  shows  the  same 

density  profiles  with  the  magnetic  field  (actually  rB  The  sharp  gradient 

0 

in  rB  is  due  to  the  fact  that  the  hot  core  increases  the  conductivity 
9  3/2 

(i.e.  o-T  )  and  prevents  diffusion  of  the  field.  It  can  also  be  seen 

Q 

tnat  the  total  current  ,  I-r(rB  ),  has  decreased  dramatically  at  the  later 

C.  0 

time  as  dictated  by  the  current  waveform. 


Clearly,  a  very  good  fit  is  obtained  for  an  ion  density  of  4x10' ^  cm-^ 
and  assumed  neon  plasma  temperature  of  165  eV.  The  other  three  calculated 
outputs  clearly  demonstrate  the  sensitivity  of  these  conditions,  relative 
to  the  line  power  outputs  produced. 

A  pinched  ion  density  of  4x10^  cm-^  with  a  diameter  of  1.8  mm  implies 
a  total  K-Shell  radiating  mass  is  34  micrograms  (yg)  per  centimeter.  The 
nozzle  puff  calculations  carried  out  by  NRL's  Plasma  Technology  3ranch 
suggest  total  puffed  masses  of  35-50  yg/cm  for  this  shot.  Therefore,  this 
preliminary  analysis  suggests  that  Gamble-II  is  very  efficient  at 
incor porating  most  of  the  puffed  mass  into  a  K-shell-radi ating  pinch.  This 
is  consistent  witn  the  very  low  L-shell  yields  that  have  been  observed  on 
this  shot  series.  It  also  contrasts  with  the  previous  BLACKJACK  5  and 
PITHON  series  where,  typically,  10-251  of  the  cold  loaded  mass  ultimately 
radiated  K-shell  photons.  This  very  interesting  and  premising  preliminary 
result  must  be  further  validated  by  more  shot  analyses  before  it  can  be 
accepted  as  fact. 
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II.  GAMBLE-II  NEON  SPECTRUM  ANALYSIS 


At  this  writing,  one  Gamble- 1 1  neon  gas  puff  shot  (if  2533)  spectrum  has 
been  reduced  to  absolute  intensities  for  the  principal  K-shell  radiating 
lines.  We  have  analyzed  this  spectrum  using  the  collisional-radiati ve 
model  described  above. 

In  shot  2533,  a  total  neon  K-shell  yield  of  2.1  kJ  was  produced.  Of 
this  2.1  kJ,  680  J  appeared  in  the  Ne  IX  1s2-1s2p1P  ("He  a")  line,  and  1040 
J  appeared  in  the  Ne  X  1s-2p  line  ("Ly  at").  Since  these  lines  obviously 
dominate  the  yield,  we  have  concentrated  our  analysis  on  reproducing  this 
power  output. 

In  each  case,  the  pinch  is  assumed  to  have  a  length  of  4  cm.  A  series 
of  four  spectra  were  calculated,  for  various  assumed  temperatures  and 
densities.  The  diameter  was  set  equal  to  1.8  mm,  in  accordance  with  the 
pinhole  imaging,  for  all  five  of  these  spectral  calculations  which 
attempted  to  fit  the  measured  power  output.  The  results  are  summarized  in 
Table  II  below. 

Table  II:  Gamble-II  Shot  2533,  Calculations  and  Measurements 
NI  Tg  He  q  power  (J/nsec)  Ly  a  power  (J/nsec) 
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3.  Same  as  Fig.  2,  except  the  quantities  are  plotted  against 
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cylindrical  plasma  diameter  for  a  sodium  ion  density  of  10  cm  at  a 
temperature  of  400  eV. 
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Fig.  2.  The  spatially  averaged  single-flight  escape  probability  P  ,  the 
collisonal  photon  quenching  probability  per  scattering  Pq,  and  the 
approximate  analytic  ultimate  escape  probability  are  plotted  against 
ion  density  for  a  400  eV  cylindrical  sodium  plasma  of  diameter  1.4 
mm.  Also  shown  is  the  exact  computed  numerical  ultimate  escape 
probability  P  . 
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50?  mass  line  where  little  or  no  disturbance  is  seen.  Thus,  a  large 

fraction  of  the  shell  (e.g.,  greater  than  or  equal  to  50?)  i3  seen  to 
remain  intact. 

The  growth  of  the  backside  of  the  shell  is  illustrated  in  Fig.  10. 

The  growth  rate  can  be  calculated  by  measuring  the  difference  between  the 

90?  mass  lines  at  Z  =  0  cm  and  at  Z  =  0.5  cm.  Both  of  these  points  are 
modes  which  are  constrained  to  move  only  in  the  R  direction.  The 
approximately  linear  growth  of  the  amplitude  is  apparent  in  this  figure.  At 
the  95?  mass  line,  the  growth  of  this  difference  is  nearly  identical  to 
that  of  the  90?  mass  line.  However,  at  the  30?  mass  line  the  growth  is 
reduced  as  shown  in  Fig.  11.  This  is  also  apparent  in  Figs.  6  and  8.  The 
growth  at  the  50?  mass  line  is  so  small  that  it  cannot  be  accurately 
measured  from  a  plot  such  as  this. 

The  growth  rate  results  are  summarized  in  Fig.  12  for  several 

wavelengths.  As  discusssed  above,  the  growth  rates  are  calculated  by 
comparing  the  radial  position  of  the  90?  mass  line  at  Z  =  0  and  Z  =  A/2. 

A  second  set  of  data  is  shown  in  which  the  procedure  for  running  the  2-D 
calculations  outlined  above  was  duplicated  except  that  the  atomic  physics 
packages  were  turned  off  in  both  the  1-D  initializing  run  and  in  the  2-D 
perturbation  runs.  These  runs  were  made  for  comparison  purposes.  The 
error  bar  denotes  a  typical  variance  in  measuring  the  growth  rate  using 
data  from  different  time  sequences. 

The  results  show  that  the  growth  rate  peaks  near  1  cm.  This  is 
approximately  the  thickness  of  the  imploding  shell.  At  long  wavelength, 
the  growth  rate  drops  off  nearly  as  the  square  root  of  A  1 .  This  is  the 
dependence  of  the  classical  Rayleigh-Taylor  growth  rate  for  a  sharp 
boundary  between  two  fluid  layers  of  semi-infinite  extent.  At  longer 
wavelengths,  the  shell  thickness  and,  consequently,  the  gradient  scale 
lengths  are  small  compared  to  the  wavelength  so  that  this  correspondence  is 
reasonable.  The  decreasing  growth  rate  at  short  wavelengths  is  also  a 
phenomena  which  has  been  seen  before  in  calculations  of  Rayleigh-Taylor 
instabilities  in  shells.1  This  can  be  attributed  to  the  fact  that  the 
instability  samples  only  a  small  portion  of  the  shell  as  contrasted  to  the 
classical  case  of  a  sharp  interface  where  the  instability  "feels"  the 
entire  density  gradient. 


49 


The  two  seta  of  data  points  shown  in  Fig.  12  exhibit  the  same 
qualitative  behavior  except  that  the  peak  growth  rate  is  lower  and  occurs 
at  a  slightly  longer  wavelength  when  the  atomic  physics  is  turned  off. 
3oth  of  these  effects  can  be  attributed  to  the  fact  that  without  the  atomic 
physics  the  thermal  energy  of  the  shell  and  background  gas  are  higher. 
Consequently ,  the  pressure  is  also  higher.  This  limits  the  compression  of 
the  shell  as  well  as  the  acceleration  of  the  shell.  The  thicker  shell 
accounts  for  the  shift  in  the  peak  growth  rate  and  the  lower  acceleration 
accounts  for  lowered  growth  rates  since  Y^„  -  /  ka  where  Y  is  the 
classical  growth  rate,  k  ;s  the  wavenumber,  and  a  is  the  acceleration.  A 
peak  growth  rate  of  1.5  x  IQ'7  sec-1  implies  a  growth  period  or  e-folding 
time  of  67  nsec.  Thus,  over  the  approximately  300  nsec  where  one  would 
expect  the  instability  to  be  active  (i.e.,  the  time  frame  over  which  the 
shell  is  accelerated  inward),  there  would  be  9.5  growth  periods.  This 
corresponds  to  an  amplitude  increse  of  90.  An  initial  amplitude  of  0.01  cm 
would  then  grow  to  0.9  cm.  This  growth  is  dramatically  reduced  at  shorter 
wavelengths  (i.e.,  high-k)  modes.  Also,  the  shell  is  unstable  over  only  a 
small  region.  Most  of  the  shell  is  stabilized  due  to  the  snowplow  effects 
of  the  magnetic  field.  The  largest  fraction  of  the  K-shell  radiation  can 
be  expected  to  come  from  the  dense  inner  part  of  the  shell  which  is  raised 
to  the  highest  temperatures  during  the  implosion.  Therefore,  these  results 
lead  one  to  expect  that  there  is  minimum  degradation  in  the  K-line  output 
from  this  effect  since  it  appears  that  at  least  50$  or  more  of  the  shell 
remains  intact. 
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Profiles  at  30-0  nanoseconds.  Peak  temperature  exceeds  300  eV. 
Shocks  propagating  radially  inward  in  core  plasma  and  radially 
outward  in  puff_?as. 
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'ig.  4.  Profiles  at  23.3  nanoseconds.  Stagnation  phase:  contact  with 

central  core  plasma  has  commenced.  Temperature  at  interface 
exceeds  100  eV. 
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Profiles  at  10.0  nanoseconds.  Coasting  phase:  forward  edge  of 
puff-gas  has  been  heated  to  a  few  tens  of  eV  through  accretion  cf 


background  plasma. 
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displacement 


Profiles  at  3.75  picoseconds.  Radial  displacement  is  given  in 
cm.,  log-jo  density  in  g/cm^,  log^Q  temperature  in  eV,  log10 

3  -10 

pressure  in  ergs/cmJ  X  10  ,  log-|Q  radial  velocity  in  cm/nsec  X 

—  •  o  - 1 4 

10  ^  and  log^  ^racl  in  ergs/cnH-nsec  X  10 


which  was  compared  to  the  self-consistent  simulation,  neglected  radiation 
but  included  ionization  dynamics.  At  certain  times  and  regions,  the  local 
energy  of  ionization  and  excitation  substantially  exceeded  the  local 
thermal  energies.  It  is  probable  that  a  pure  ’’hydrodynamic"  calculation 
with  simplified  prescriptions  for  ionization  energy  and  average  ionization 
state  would  produce  a  substantially  different  plasma  evolution. 

Since  radiation  plays  such  a  crucial  role,  it  must  be  generated  and 
transported  accurately.  As  we  have  shown,  optically  thin  radiation  ignores 
important  transport  of  energy  by  radiation  and  grossly  overestimates 
losses.  LTE  or  coronal  models  for  the  level  populations  will  also  lead  to 
inaccurate  energetics.  We  found  that  the  atomic  populations  in  the  regions 
of  pea*  radiative  emission  were  substantially  different  from  corresponding 
LTE  or  coronal  populations,  in  some  cases  by  orders  of  magnitude. 

The  neon  model  which  we  employed  was  a  relatively  modest  one,  with 
only  27  atomic  levels  and  13  emission  lines.  We  have  available  more 
extensive  neon  models  with  nearly  a  hundred  levels  and  with  hundreds  of 
lines.  How  much  accuracy  has  been  given  up  by  using  a  reduced  model? 
Experience  has  shown  that  simple  models  can  underestimate  radiative  losses 
bv  perhaps  a  factor  of  two.  Thus,  we  plan  to  perform  an  identical 
calculation  with  a  full  neon  model,  and  it  is  expected  that  the  radiation 
yield  will  be  even  higher. 
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Corresponding  profiles  of  density,  temperature,  pressure  and  velocity 
are  shown  in  Figs.  13  and  14.  At  30.0  nanoseconds  (Fig.  13),  the  profiles 
are  very  similar  in  the  interface  region.  However,  the  role  of  radiation 
in  transporting  energy  is  clearly  evident.  The  dense  core  is  still  very 
cold  compared  with  the  complete  simulation  with  radiation  transport.  Also, 
the  temperatures  through  the  forward  half-centimeter  of  the  puff-gas  plasma 
are  a  factor  of  about  3  lower  than  in  the  full  simulation.  A  comparison 
with  the  profiles  in  Figs.  5  through  7  suggest  how  greatly  "thin"  radiation 
can  overestimate  energy  loss. 

In  spite  of  the  differences  in  the  calculations,  assembly  occurs  at 
about  the  same  time.  The  situation  at  50.0  nanoseconds  is  shown  in  Fig. 

14.  There  are  qualitative  similarities  with  Fig.  7,  but  substantial 
differences  are  also  evident.  Temperatures  and  densities  differ  by  at 
least  a  factor  of  two  almost  everywhere.  Also,  a  substantial  portion  of 
the  plasma  is  moving  radially  outward  by  this  time. 

A  corresponding  energy  history  of  this  simulation  is  shown  in  Fig. 

15.  The  early  history  is  similar  to  that  shown  in  Fig.  10.  3y  30 

nanoseconds,  however,  it  is  evident  that  substantially  more  energy  resides 
in  thermal  and  ionization  energies.  The  relative  minimum  in  kinetic  energy 
occurs  earlier  and  is  not  as  sharp.  As  can  be  seen  from  the  radiation 

curve,  if  radiation  post-processing  had  been  employed,  all  of  the  available 
plasma  energy  would  have  been  lost  within  the  first  30  nanoseconds. 

IV.  CONCLUSIONS 

A  puff-gas  imploding  onto  a  central  core  plasma  has  promise  as  a 

"eHable  way  to  achieve  sufficiently  high  density  and  temperature  over 
sufficiently  large  volumes  and  timescales  to  efficiently  convert  kinetic 
energy  to  radiation.  In  the  case  we  have  studied,  a  neon  puff-gas 
implading  ante  a  neon  core  plasma,  approximately  two-thirds  of  the  total 
plasma  energy  was  converted  to  K  and  L-shell  radiation  by  100  nanoseconds. 

3e cause  the  radiative  energy  loss  and  the  transport  of  radiation  in 

the  plasma  substantially  modify  the  hydrodynamic  response,  a  self- 

consistent  calculation,  linking  radiation  and  fluid  transport  with  the 
ionization  iynamics,  must  be  performed.  The  "hydrodynamic"  calculation, 
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The  energy  history  of  the  plasma  throughout  the  first  120  nanoseconds 
is  shown  in  Fig.  10.  The  energy  in  various  categories  is  shown  as  a 

C  11 

function  of  time,  on  a  logarithmic  scale  from  10J  to  10  ergs/cm. 

Initially,  about  98J  of  the  total  energy  is  kinetic;  the  other  2%  is 
in  thermal  and  ionization  energies.  The  energy  radiated  away  begins  to 
increase  sharply  just  before  20  nanoseconds,  as  the  puff  gas  hits  the 
central  core.  Thermal  and  ionization  energies  also  begin  to  increase 
sharply  at  this  time.  Kinetic  energy  goes  through  a  well  defined  minimum 
at  about  60  nanoseconds.  Althrough  peak  compression  on  axis  occurs  at 
about  50  nanoseconds,  substantial  portions  of  the  puff  gas  plasma  continue 
to  move  inwards  until  about  60  nanoseconds.  Notice  that  by  60  nanoseconds, 
almost  half  of  the  total  plasma  energy  has  been  radiated  away.  3y  200 
nanoseconds,  about  2/3  of  the  total  energy  has  been  converted  to 
radiation.  This  makes  it  very  difficult  to  justify  a  "hydro-only" 
calculation,  or  a  calculation  where  the  radiation  energetics  are  not 
carefully  performed. 

Figure  11  shows  the  calculated  emission  spectrum  from  the  plasma  at 
23.3  nanoseconds,  just  after  the  puff-gas  has  contacted  the  central  core. 
The  energy  range  is  10  eV  to  10  keV,  and  the  spectral  intensity  is  in  units 
of  ergs/cm-keV-sec.  At  this  time,  there  is  substantial  L-shell  emission, 
but  little  K-shell  radiation.  Most  of  the  L-shell  radiation  is  coming  from 
the  denser,  cooler  portion  of  the  interface  region,  while  the  K-shell 
radiation  is  coming  mostly  from  the  hotter  but  less  dense  portion  just 
outside  the  core  plasma. 

At  50.0  nanoseconds  (Fig.  12),  the  emission  spectrum  shows 
substantially  more  power  in  the  K-shell  and  the  L-shell.  A  broader  annulus 
of  plasma,  situated  outside  the  core  plasma,  is  producing  most  of  the 
radiation.  Although  the  very  dense  plasma  on  axis  is  radiating  strongly, 
its  volume  is  small,  and  the  temperature  is  relatively  low. 

A  similar  calculation  was  performed  without  radiation  cooling,  but 
with  detailed  atomic  physics,  so  that  the  ionization  energies  and  average 
ionization  numbers  would  be  sel f-consistently  calculated.  Although  the 
radiative  cooling  was  turned  off,  optically  thin  radiation  was  calculated 
for  comparison  with  the  simulation  discussed  above. 


inwards.  Some  radial  compression  of  the  puff-gas  is  in  evidence.  Thermal 
conduction  has  produced  a  nearly  isothermal  region  between  the  outer  edge 
of  the  core  plasma  and  the  forward  millimeter  of  the  puff-gas.  Radiative 
effects  at  this  time  are  negligible. 

At  t  =  23-3  nanoseconds  (Fig.  U) ,  bhe  puff-gas  is  in  full  contact  with 
the  core.  The  compression  at  the  interface  has  created  a  large 
over pressure ,  and  the  temperature  in  this  region  exceeds  TOO  eV.  A  narrow 
region  of  radiative  cooling  is  situated  exactly  at  the  interface.  Adjacent 
regions  are  net  absorbers  of  radiation  -  they  are  being  heated  by  the 
absorption  of  soft  x-rays  produced  near  the  interface. 

At  t  =  30.0  nanoseconds  (Fig.  5),  the  peak  temperature  near  the 
interface  exceeds  300  eV.  The  large  overpressure  has  caused  shocks  to 
propagate  radially  inward  in  the  core  plasma,  and  radially  outward  in  the 
puff-gas.  Heating  from  the  soft  x-rays  produced  in  the  interface  region  is 
evident  in  the  dense  core.  Radiation  heating  is  also  taking  place  deep  in 
the  puff  gas.  The  effects  of  radiation  transport  will  be  more  clearly  seen 
when  comparison  is  made  with  hydro-only  profiles. 

At  -40.0  nanoseconds  (Fig.  6),  the  shocks  are  more  fully  developed. 
Thermal  conduction  is  creating  a  nearly  isothermal  region  near  the  leading 
edge  of  the  puff-gas. 

At  50.0  nanoseconds  (Fig.  7),  the  inward-going  shock  has  reached  the 
origin.  Although  the  density  approaches  solid  density  on  axis,  the  plasma 
temperature  there  is  still  only  a  few  tens  of  eV. 

By  60.0  nanoseconds  (Fig.  8),  portions  of  the  plasma,  including  the 
dense  core,  are  moving  "adially  outwards.  The  outward  going  shock  is  still 
in  evidence,  and  the  outer  edge  of  the  puff-gas  is  near  its  minimum  radius 
( a  bo  ut  0.5cm). 

Finally  70.0  nanoseconds  (Fig.  9),  the  entire  plasma  is  moving 
radially  outwards.  The  core  plasma  remains  cool,  3ince  the  density  is 
still  high  enough  to  retard  thermal  conduction.  Furthermore,  the  radiative 
emission  of  this  regio-  is  in  approximate  balance  with  the  radiative 


absorot ion . 


where  Fp^  is  the  rate  of  energy  loss  in  zone  k  due  to  a  discrete  radiative 
process  (or  frequency  group)  P,  and  Cp(<j  is  the  radiative  coupling  of  zone 
k  to  zone  j  for  that  process.  The  couplings  are  functions  of  opacity, 
integated  over  process  and  photon  path.  In  the  probabilistic  model,  a 
matrix  of  couplings  must  be  computed  for  each  bound-bound,  bound-free  and 
free-free  process;  for  the  multifrequency  model,  they  must  be  computed  for 
each  discrete  frequency.  In  this  way,  the  net  cooling  and  heating  by 
radiation  emission  and  absorption  among  the  various  zones  of  the  plasma  is 
accounted  for  accurately. 


III.  RESULTS 


We  n^ve  performed  a  calculation  involving  a  cylindrical  neon  puff  gas 

of  density  5x10-°g/cm^,  with  inner  radius  0.55  cm  and  outer  radius  1.95 

cm,  imploding  radially  at  a  velocity  of  3x10^cm/sec  on  a  uniform  neon 

central  core  plasma  of  density  5x10~3g/cm3  anq  ra(jiUS  o.10  cm.  The  initial 

configuration  is  as  shown  in  Fig.  1.  The  temperature  of  the  puff  gas  is 

taken  to  be  5 eV  initially,  and  the  core  plasma  temperature  is  about 

.O^eV.  There  is  assumed  to  be  a  tenuous  background  plasma  in  between  the 
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puff  gas  and  core  plasmas,  of  density  5x10  g/cm-5. 

Profiles  of  density,  temperature,  total  pressure,  velocity  and 
radiative  cooling  at  various  times  in  the  calculation  are  given  in  Figs.  2 
through  9.  The  ordinate  in  all  of  the  plots  is  logarithmic,  ranging  from 

_  Q  +.  ii 

10  to  10  .  The  various  plasma  profiles  are  suitably  scaled  so  that  they 

can  be  plotted  together:  density  in  g/cm^,  temperature  in  eV.,  pressure  in 

■3-12  -2 

ergs/cm-1  x  10  ,  velocity  in  cm/nsec  x  10  and  radiative  cooling  in 

er gs/ cm J-n3ec  x  10  .  The  abscissa  is  linear,  and  ranges  from  0.0  (the 

central  axis)  to  1.0  cm. 

At  t  *  3-75  picoseconds  (Fig.  2),  the  plasma  is  very  close  to  the 
initial  configuration.  There  is  a  slight  warming  of  the  plasma  evident  at 
the  background-puff  gas  interface.  Note  that  only  the  inner  1.0  cm  of  the 
calculation  is  shown;  the  puff  gas  initially  extends  to  1.95  cm. 

At  t  =  10.0  nanoseconds  (Fig.  3),  the  forward  edge  of  the  puff-gas  has 
been  heated  to  a  few  tens  of  eV  through  accretion  of  the  background 


plasma.  At  this  time,  the  puff-gas  plasma  is  essentially  coasting  radially 


c.  Radiation  Emission  and  Transport 


Radiation  emission  from  a  plasma  and  its  opacity  are  dependent  on  the 
local  atomic  level  population  densities.  Except  for  optically  thin 

Plasmas,  however,  the  level  populations  depend  on  the  radiation  field, 
since  optical  pumping  via  photoio.ni zation  and  photoexcitation  can  produce 
significant  population  redistr i bution .  Thus,  the  ionization  and  radiation 
transport  processes  are  strongly  coupled  and  must  be  solved  self- 
consistently .  In  this  model,  an  iterative  procedure11  is  used,  where  level 
populations  are  calculated  using  the  radiation  field  from  the  previous 
iteration,  then  using  these  populations  to  calculate  a  new  radiation  field 
and  recalculating  populations  until  convergence  is  reached. 

Three  distinct  radiation  transport  schemes  have  been  developed,  and 

can  be  used  interchangeably  in  the  code:  a  probabilistic  scheme,  a 

mul ti frequency  scheme  and  a  hybrid  method.  Descriptions  of  these  transport 
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methods  have  appeared  elsewhere.  ’  J  The  probabilistic  model  was  employed 
in  the  present  calculation,  since  the  total  plasma  mass  is  relatively 
small . 

1  2 

The  probabilistic  model  forms  local  angle  and  frequency  averaged 
escape  probabilities  for  each  emission  line  and  for  each  bound-free 
process .  Free-free  radiation  is  treated  with  a  multi  frequency  formalism. 
The  radiation  transport  and  emission  spectra  are  calculated  from  these 
escape  probabilities.  The  method  can  treat  comprehensive  atomic  models  and 
provides  good  overall  energetics,  but  cannot  calculate  certain  spectral 
details  and  becomes  inaccurate  at  very  high  optical  depths. 

In  all  of  the  models  outlined  above,  local  ionization  state-dependent 

inner-shell  opacities  are  included,  since  these  processes  are  very 

important  in  the  cool,  dense  plasma  regions.  Inner-shell  photoionization 

ci-'oss  sections  for  the  neutral  element  are  taken  from  the  fits  by  3iggs  and 

lighthili,14  and  the  positions  of  the  ionization-dependent  absorption  edges 

1  s 

are  taken  from  the  Hartree-Fock  calculations  of  dementi  and  Roetti.  J 

The  local  rate  of  energy  change  in  zone  j,  due  to  radiation  transport 
is  given  by 


t  . 

J 


(9) 
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b. 


Ionization  and  Atomic  Physics 


The  ionic  populations  in  the  plasma  may  be  characterized  by  a  set  of 
atomic  rate  equations  of  the  form 


df . 
1 

dt 


E  W . .  f. 

j  J1  J 


E  W.  .  f. 
i  1J  1 


(6) 


where  f^  is  the  fractional  population  of  atomic  level  i,  and  W j ^  is  the  net 
reaction  rate  describing  the  transition  from  initial  state  j  to  final  state 
i.  An  equation  of  this  type  is  constructed  for  each  of  the  atomic  levels 
included  in  the  model. 


For  sufficiently  dense  plasmas  of  the  sort  we  expect  to  model  with  the 
implosion  code,  the  effective  populating  and  depopulating  rates  are 
generally  fast  compared  with  the  hydrodynamic  response.  An  equilibrium 
assumption  can  be  justified,  which  involves  dropping  the  explicit  time 
dependence  in  Eq .  (6).  The  plasma  is  then  said  to  be  in  collisional- 

radiative  equilibrium  (CRE),^  whereby  the  plasma  ionization  state  responds 
instantaneously  to  changes  in  hydrodynamic  quantities. 

The  rate  coefficients  that  are  used  to  calculate  the  populating  and 

depopulating  rates,  Wj^,  are  calculated  using  various  scattering 

techniques.  The  processes  included  in  this  calculation  and  the  methods 

used  in  calculating  the  corresponding  rate  coefficients  are  summarized 
Ii_  i  n 

elsewhere. 

Once  the  set  of  rate  equations  (including  the  radiation  transport)  has 
been  solved  for  the  level  populations,  f^,  the  electron  density  can  be 
cal culated , 

N  =  E  z.  f.  N,  (7) 

e 


where  z-  is  the  ionic  charge  of  level  i  and  Nj  is  the  total  ion  density. 
The  ionization  and  excitation  energy  can  also  be  calculated  by 


e . 
i 


f . 

l 


(3) 


wnere  y.  is  the  energy  of  level  i,  measured  from  the  ground  state  of  the 
neutral  atom. 
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wr.ere  p  is  mass  density,  u  is  velocity,  P  is  pressure,  is  total  energy 
density,  £„ad  is  the  rate  of  energy  loss  or  gain  due  to  radiation,  n  is  the 
tnermal  conductivity,  and  N  is  the  ion  density.  The  thermal  conduction  is 
calculated  implicitly,  using  an  iterative  Crank-Micholson  scneme. 

Since  iensities  did  not  exceed  solid  density  in  this  study,  a  simple 
equation  of  state  was  assumed  with 


D 


1 

5 


o 

pu“  - 


<d) 


where  is  the  potential  energy  due  to  ionization  and  excitation.  A  non¬ 
ideal  equation  of  state  taking  account  of  ionization  energy  and  degeneracy 
pressure  is  employed  in  cases  where  the  density  exceeds  solid  density.  A 
single  temperature  model  was  employed, 

kT  =  - - - —  ,  (5) 

(p/m  )  (1+Z) 

where  m-r  is  ion  mass,  and  T  is  temperature.  The  ionization  energy,  ,  and 
effective  charge,  Z  are  calculated  from  the  ionization-radiation  equations 
which  are  explained  below.  A  single  temperature  assumption  is  valid  in  the 
core  plasma,  where  the  equilibration  time  is  of  the  order  of  picoseconds, 
and  it  is  adequate  in  the  stagnation  region,  where  the  equilibration  time 
can  be  of  the  order  of  nanoseconds.  In  the  blowoff  plasma,  it  is  a 
marginal  approximation,  but  the  consequences  are  small,  since  little 
radiation  is  emitted  from  this  region,  and  since  most  of  the  thermal  energy 
is  carried  by  the  electrons  in  the  blowoff. 

The  local  rate  of  change  of  energy  due  to  radiation  transport,  e  , , 

"ad 


will  be  discussed  below. 


Implosion  of  a  Neon  Puff-Gas  Plasma  Onto  a  Low-Density  Plasma 


I.  INTRODUCTION 


Z-pineh  implosions  of  puff-gas  plasmas  onto  low  density  foams  are 
presently  being  carried  out  at  Sandia  National  Laboratory.  The  collapse  of 
an  imploding  puff-gas  on  a  central  core  produces  hot  moderate  density 
plasma  over  reasonably  large  volumes  and  timescales.  Sandia  reports  high 
"adiative  conversion  efficiencies  from  relatively  clean,  reproducible 
implosions.  We  have  simulated  a  typical  Sandir  Z-pinch  experiment  with  a 
self-consistent  ID  radiation-hydrodynamic  scheme.  Such  simulations  can 
help  in  understanding  the  underlying  physics  of  these  experiments,  and  can 
lead  to  diagnostic  and  predictive  capability. 


II.  THEORETICAL  MODEL 


Imploding  a  cylindrical  puff-gas  plasma  onto  a  central  core  can  result 
in  the  radiative  loss  of  a  substantial  fraction  of  the  total  plasma 

energy.  The  radiative  emission  depends  strongly  on  the  local  plasma 

density  and  temperature.  In  addition,  optical  pumping  and  energy  transport 
by  photons  of  all  frequencies  can  influence  the  plasma  temperature  and 
degree  of  ionization,  and  can  modify  the  hydrodynamic  evolution  of  the 
plasma.  Thus,  the  hydrodynamic  development  and  atomic  physics  of  the  puff 
gas  and  core  plasmas,  as  well  as  the  transport  of  radiation,  must  be 
calculated  self-consistently . 

Discussion  of  the  theoretical  model  can  be  separated  for  convenience 
as  follows:  (a)  hydrodynamics  and  thermal  conduction  (b)  ionization  and 

atomic  physics  and  (c)  radiation  emission  and  transport. 

a.  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  mass,  momentum,  and  total  energy 

are  transported  in  one  dimension  using  a  numerical  scheme  with  a  sliding- 

2 

zone  version  of  f lux-corrected  transport.  A  special  gridding  algorithm  is 
used  which  moves  zones  in  a  Lagrangian  fashion  and  adjusts  the  mesh  in 
order  to  resolve  steep  gradients  in  the  flow.  The  hydrodynamic  equations 
solved  are 


hydro  or ofiios 
40.  CO  NSEC. 
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Fig.  6.  Profiles  at  ^0.0  nanoseconds.  Core  plasma  heated  by  soft  x-rays 
from  interface  region.  Radiation  energy  transport  evident  in 
puff-gas  out  to  about  O.d  cm. 
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displacement 


Profiles  at  50.0  nanoseconds.  Assembly  phase:  inward 

propagating  shock  arrives  at  axis;  peak  density  approaches  solid 


tumpGruLure 


gas  plasmas  moving  radially  outward.  Global  kinetic  energy  near 


•I.M  •I.ODt  -1.77*  LW 

(«•  M  l-f  l*U  »t> 

IW*  LUflM-M 

Emission  spectrum  at  23-3  nanoseconds.  Spectral  intensity  (from 
10'5  to  lO^1  ergs/sec-cm-keV )  is  plotted  as  a  function  of  photon 
energy  (f’-'om  10  eV  to  10  keV).  At  this  time,  L-shell  emission 

dominates. 
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Fig.  13.  Profiles  at  30.0  nanoseconds  neglecting  effects  of  radiation. 

The  core  plasma  remains  very  cold,  and  there  is  no  evidence  of 
heating  deep  in  the  puff-gas. 
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14.  Profiles  at  50.0  nanoseconds  neglecting  effects  of  radiation. 
Assembly  on  axis  occurs  at  about  the  same  time,  but  there  are 
substantial  differences  in  the  plasma  density  and  temperature 
profiles. 


energy  history 


thin  rca 


0.  0  48.  0  96.  0 


time  (nsec.  ) 

Fig.  15.  Energy  history  of  puff  gas  implosion  neglecting  effects  of 
radiation.  Optically  thin  radiation  is  calculated  but  it  does 
not  remove  energy  from  the  plasma. 
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E .  Velocity  Scaling  of  the  Radiatively  Dominated  Z-Pinch 


I.  THE  ELECTRODIFFUSIYE  MODEL 

For  a  simple  diode  the  general  formulation  of  a  fast  Z-pinch  implosion 
model1  can  be  schematically  represented  as  shown  in  Fig.  1.  The  input 
current  waveform  at  some  radius  and  the  initial  plasma  load  conditions,  in 
any  of  several  possible  representations ,  provide  a  unique  initial  value 
problem.  The  complete  electromagetic  calculation  being  quite  difficult,  a 
common  approximation  is  that  of  neglecting  all  displacement  currents  and 
the  details  of  capacitive  energy  storage. 

a.  Diffusive  Approximation 

If  one  gives  up  some  information  concerning  the  details  of  the  diode 
fields  and  makes  the  assumption  that  the  incoming  and  outgoing  wave 
components  are  in  detailed  balance,  then  the  Hertz  wave  equation  (for 
Z  (r,t))  can  be  transformed  to  a  diffusion  equation^ 

~  E  -  {u_13  u3  E}  7-!-=  -  E3  Sin  l  -  3  E..  -  B3  0,  (1) 

Dx  u  u  4ttE  x  l  t  th  x 


where  Eth  is  the  thermoelectric  field,  L  is  a  dimensionless  con¬ 
ductivity,  0  -  V  /c,  E  -  E  ♦  0B,  E  =»  E  +  0B  ♦  E  ,  D/Dx  »  3  ♦  03  , 

I  i.  Li  X  Q  in  X  u 

and  u  =  r/r  ,  x  -  t/t  ,  .  -  ct  .  The  fluid  resDonse  to  the  electro- 

o  o  o  o 

magnetic  stresses  and  heating  is  embodied  in  the  relations 
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In  these  expressions  9  T  +  (2/3]eT/z,  T  =  Tt  +  zT  ,  in  =  Mr  +  z  m  ; 

e  e  I  I  e  I  max  e 

X  ,  ,  is  the  thermal  conductivity,  x  the  plasma  relaxation  time,  E  the 
0  X  0  * 

ambipolar  radial  field  (with  p  its  induced  charge  density),  and  QraC[  is 
the  net  (local)  radiative  heating  or  cooling.  The  radial  electric  field  is 
established  as  a  solution  to  an  integral  equation  derived  from  the  radial 
component  of  Ohm’s  law.  Any  drift  speed  limitation  on  current  i3  supple¬ 
mented  by  a  (nonlinear)  change  in  Z  where  the  local  E  field  requires  it. 

b.  Discrete  Realization 


Notation :  The  (nonuniform)  time  levels  for  any  variable  are  indexed 

with  a  leading  superscript;  the  Lagrangian  fluid  mesh  is  denoted  by  r,  and 

its  material  derivatives  by  a  superscripted  dot  (or  dots).  Spatial 

indexing  is  denoted  by  trailing  subscript. 

Method ;  The  base  representation  for  all  the  hydrodynamic  variables  is 

the  set  (  N.,  LT  ,  x9  .,  lr .  [  and  j1  E  .f,  the  diffusing  electric 

c  j  c  Ij  c  ej  j  1  c  z , j  ’ 

field.  The  relationships  among  the  basic  fluid  variables  on  the  Lagrangian 

mesh  are  illustrated  in  Fig.  2.  The  numbers  of  cell  ions  ;  N.}  are  a 

2  J 

conserved  vector  of  ions/cm  resident  in  the  (compressible)  cell  and 
assigned  spatial  locations  given  by  the  cell  center  position,  the  equal 
area  point. 

In  general,  at  any  subcycle  point  t*  within  ( ^ t ,  1  V],  one  must  map 

*  ,  |  *  .  *.  *  .  , 

Q  .  J  -►  1  Tt.  ,  E  .  ,  Q  .  in  such  a  way  that  the 
c  ej  ‘c  Ij  c  zj  c  ej 

. »  *  *  *  . 

intermediate  fields'  1  r  j ,  V  .,  S  ,  J  .  are  all  consistent  witn  one 

rj  c  9k’  zj  ,  * 

another.  Only  the  time  can  enter  as  a  generator  of  r^  and  vr j >  because 


t  ,  Tt. ,  E  . , 
'  c  Ij  c  zj 


89 


or'  tne  Lagrangian  nature  of  the  mesh. 

*  * 

Although  the  V  .  and  3^.,  calculations  must  be  considered  coupled 
rj  c  0j 

problems,  a  magnetic  field  calculation  from  any  given  velocity  profile  is 
the  pivotal  element  in  achieving  the  desired  self-consistency. 

If  tne  mesh  and  its  velocity  are  presumed  given,  or  at  least 
provisionally  established  subject  to  iterative  refinement,  then  there  is 
just  one  current  density  profile 
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with  =  (2/cr)  (2trr  1  J^)dr1  ,  and  x  =  ~  (1.7588  •  IO^t  )B  ,  a  nonlinear 

integral  equation,  using  Bragi.nsk i i ' s  functions  8,5.  This  information  is 

.  .  * 

sufficient  to  define  a  derivative  function,  Y.  =  Y.(Y,  ,t  )  with  Y,  = 

♦  4  «  J  J  *  K 

;  r  ,  S  ,  3  ,  ;  and  the  Jacobian  is  either  lower  Hesse.nberg  (no 

IK  ^  ZK  0  0  K 

photons)  or  a  full  matrix  (with  radiation  transport). 

The  mesh  t.  ,  rv  and  Y^  are  then  integrated  by  a  fully  implicit  NML 
technique-1,  using  an  estimated  integrating  factor  to  track  the  "stiff" 
elements  of  the  set.  Each  location  Yj  thus  obtains  a  local  filter  for 
integration  which,  for  a  2nd  order  method,  is 
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b.  Limitations  and  Remedies 

Once  the  displacement  current  has  been  neglected  one  must  also  approx¬ 
imate  in  some  way  the  source  of  field  momentum  (g  =  -S  )  and  the  proper 

r  c  r 

formulation  is  that  of  a  flux  limit.  This  flux  limit  takes  the  form  of  a 
filter  derived  from  the  detailed  balance  of  field  and  fluid  momentum 
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which  implies  a  bound  on  J  dr  O^g)  over  each  resolution  domain  for 
/  dr  OP)  as  determined  by  the  artificial  viscosity  model.  The  result  is 
a  generalized  Rankine-Hugoniot  jump  condition 


*22  - 
p»E  +  B  +n.  T  (total  pressure).  T  -  zT  +  T,  (effective  fluid 
ion  r  e  ion 

temperature).  pQ  -  Mion  (density  at  cell  center).  The  flux  limit  arises 
as  a  bound  on  the  spatial  gradient  of  3  admitted  in  the  E  field  source 


term  33..B,  i.e. 
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The  limit  prevents  unstable  growth  of  Ez  on  the  outer  surface  of  the  model 
plasma  due  to  underresolved  flow  gradients. 

A  second  and  more  fundamental  limitation  arises  due  to  acausal 
propagation  of  magnetic  disturbances.  This  may  obtain,  for  example,  when 
an  interior  zone  increases  its  conductivity,  due  to  compression  or  Ohmic 
heating,  and  draws  more  current.  The  added  current  will  produce  an 

immediate  rise  in  strength  of  B  in  all  exterior  zones,  but  the  resulting 
change  in  the  (magnetically  insulating)  Sr3g  component  of  the  total  E' 
field  is  not  a  diffusive  source  of  E  in  those  exterior  domains.  When  this 
source  does  not  act  the  flow  appears  "over-insulated",  the  total  E'  field 
(in  the  plasma  frame)  reverses  and  so  must  J  in  the  exterior  domains.  In  a 
proper  electromagnetic  calculation  this  source  would  be  reflected  in  3t<J, 
the  E  field  would  propagate  inward  rapidly  and  cancel  the  change  in  the 
BB.  The  current  reversal  would  not  be  as  likely  to  occur.  For  fixed 
3^J  the  wave  equation  source  term  produces  a  first  integral 


3tE 


-  V 


rev 


While  this  source  and  its  first  integral  are  unavailable  to  the  diffusive 

theory,  an  explicit  source  for  E  is  added  which  turns  on  when  a  reversal  of 

J  occurs  and  seeks  to  drive  the  lab  frame  E  in  a  direction  that  will  quench 

the  reversal.  In  the  above  relation  the  integration  interval  t  forms  a 

rev 

natural  free  parameter  with  which  to  calibrate  the  "battery"  source  term  to 
a  time  scale  below  that  typical  of  magnetosonic  wave  transit  and  limit  the 
resulting  stiffness  to  acceptable  levels. 

A  third  practical  limitation  to  the  model  is  the  difficulty  in 
obtaining  current  and  voltage  measurements  (at  the  outer  radius)  which  are 
well  synchronized  and  comparably  accurate.  If  the  experimental  voltage  and 
current  traces  have  any  relative  phase  error,  the  subtraction  of  these 
waveforms  to  produce  a  plasma  voltage  is  subject  to  excessive  noise  and 


tends  to  produce  rather  wild  boundary  conditions  for  the  model. 

The  most  effective  remedy  to  date  is  to  insist  on  continuity  in  the 
waveguide  impedance  at  the  plasma  boundary  and  let 
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where 


E  =  ghost  point  boundary  value. 

Ip  =  axial  current  flowing  in  the  plasma, 

I  =  known  current  waveform  near  the  load, 
cEj  -  Ez  value  at  the  outermost  cell  center. 

II.  VELOCITY  SCALING  BEHAVIOR 


a.  A  simple  diagnotic 


The  study  consists  of  driving  a  constant  initial  configuration  with  a 
current  pulse 

i  {  <t-t  >  \ 

I  (t)  -  fl  -5  (  1  ♦  tanh  ■  ■  °  )  ;  (5) 

gen  o  2  \  tg  / 

for  a  load  mass  (in  Neon)  of  63  yg/cm,  a  peak  initial  density  at  1.25  cm, 
and  a  current  IQ  typical  of  Gamble  II  currents,  the  load  implodes  in  about 
125  ns.  The  value  of  "f"  defines  the  scale  parameter. 

In  order  to  compare  the  1-D  calculation  with  the  0-D  snowplow  model\ 
one  may  mass  average  the  velocity  and  electric  fields  as 

V.  *  m(  l  N<  V2.  /  IN.)172  (6a) 

load  l  icri  “  i  J 
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and  from  E develop  the  expected  scale  velocity  arising  from  the  snowplow 
model 


where  pQ  is  the  initial  average  mass  density.  The  ratio  of  these  two 
velocities  generally  represents  a  rough  comparison  of  the  1-D  model  pinch 
speed  and  that  of  the  snowplow,  and  tends  to  peak  at  a  time  near  the  middle 
of  the  Z-pi.nch's  run  down. 


H  a)  = 
rad , 


implosion 


/  V  ,  ) 

snowplow 


and  one  may  examine  the  quality  of  the  snowplow  model  as  a  predictor  of 
pinch  speed  (and  available  kinetic  energy  for  later  radiation). 


b.  Conclusion 


As  shown  in  Fig.  3.  the  snowplow  scaling  with  f:  0.15  <  f  <  1.75  is  by 
no  means  perfect.  The  effect  of  early  radiaton  on  the  scaling  error  is 
most  acute  for  the  low  current  cases.  However,  the  general  trend  is  to 
reduce  the  pinch  velocity  and  narrow  the  load  width,  as  shown  in  the  fluid 
trajectory  comparison;  Fig.  4a, b.  At  the  higher  current  levels,  the 
closing  of  the  S(f)  and  S(f,  Qrad  -  0)  appears  to  be  asymptotic,  rather 
than  showing  a  cross  over. 

A  more  extensive  study  can  thus  indicate  the  corrections  needed  in 
V3nowpiow  to  use  ^  33  a  Sui-de  *n  predicting  the  late  time  available 
kinetic  energy  resident  in  the  load. 
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F.  Ion-beam  deposition,  heating,  and  radiation  from  an  aluminum  plasma 


I.  INTRODUCTION 

Recent  advances  in  pulsed  power  and  diode  technologies  have  made  it 
possible  to  deliver  intense  proton  beams  onto  a  planar  target  with  high 
efficiency.  Thus,  the  interaction  of  energetic  charged  particle  beams  with 
materials  of  different  composition,  density,  temperature,  ionization  state, 
etc.,  has  become  an  active  field  of  research  (for  example,  see  citations  1- 
16  and  their  included  references).  One  significant  application  of  this 
research  is  the  use  of  ion  beams  to  generate  x-rays  for  materials  research; 
other  applications  include  ion  implantation  studies  and  ICF  investigations. 

Many  studies  of  intense  ion  beam-materials  interactions  have 
introduced  a  variety  of  assumptions  and  approximations  to  characterize  the 
atomic  physics,  ionization  dynamics,  absorption/ emission  physics,  and 
radiative  transfer.  For  example,  Saha  equilibrium  is  often  invoked,  and 
radiation  transport  is  treated  via  a  photon  diffusion  aproximation  with  use 
of  Planck  (black-body)  radiation  distributions.  These  treatments  may 
prevail  in  the  cooler  denser  regions  but  not  for  the  entire  plasma.  In 
this  paper,  a  collisional  radiative  equilibrium  (CRE)  treatment  of  the 
ionization  balance,  a  detailed  configuration  accounting  (DCA)  of  the  atomic 
level  structure,  and  a  radiation  transport  algorithm  based  on  a  combination 
of  probabilistic  and  frequency-by-frequency  ray  tracing  techniques  are 
used.  The  CRE  model  generates  Saha-Boltzmann  populations  in  the  denser, 
cooler  regions  of  the  plasma  and  reduces  to  the  corona  model  in  the 
rarefied,  hotter  regions.  The  radiation  transport  algorithm  reproduces  the 
results  of  the  photon  diffusion  approximation  in  the  cooler,  denser  plasma 
regimes.  Thus,  the  approach  taken  in  this  paper  is  more  general  than  those 
which  employ  the  foregoing  assumptions,  but  it  reproduces  their  results  in 
their  regimes  of  validity;  moreover,  our  approach  is  applicable  in  regimes 
where  these  assumptions  cannot  be  used.  This  model  has  been  described  in 
detail  elsewhere.1^  Instead  of  using  the  Bethe  or  LSS  formalism1^  for 
calculating  stopping  power  by  the  bound  electrons  in  the  material,  the 
detailed  generalized  oscillator  strength  (GOS)  stopping  cross  sections  of 
McGuire,  Peek,  and  Pitchford1^  for  AS.  ions  are  used. 
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The  sequence  of  events  for  charged  particle  beam-matter  interaction 
involves  consideration  of  several  phenomena.  First,  there  is  the  energy 
deposition  within  the  slab  by  the  beam.  Then,  there  is  the  subsequent 
hydrodynamic  response  of  the  slab  as  it  responds  to  this  energy  influx;  a 
blow-off  region  develops  on  the  front  (beam)  side  of  the  plasma,  while  the 
back  surface  recoils  away  from  the  beam.  As  the  slab  plasma  heats  up,  the 
plasma  ions  emit  radiation  which  can  affect  the  energy  balance  of  the 
plasma.  Thus  at  any  instant  in  the  plasma,  there  is  the  interplay  of 
deposited  beam  energy,  thermal  conduction,  hydrodynamic  motion,  radiation 
transport,  plasma  internal  energy,  magnetic  fields,  etc.  to  be  taken  into 
account . 

With  the  exception  of  magnetic  and  electric  field  effects,  this  paper 
represents  an  attempt  to  treat  all  these  effects  self-consistently  in  a 
one-dimensional  simulation  of  a  proton  beam  incident  on  an  initially  solid 
density  A1  slab. 

II.  THEORETICAL  MODEL 

The  interaction  of  tne  ion  beam  with  the  target  is  strongly  dependent 
on  the  local  temperature  and  degree  of  ionization  in  the  plasma.  In 
addition,  optical  pumping  and  energy  transport  by  photons  of  all 

frequencies  can  influence  the  plasma  temperature  and  degree  of  ionization, 
and  can  modify  the  hydrodynamic  response  of  the  slab.  Thus,  the  deposition 
of  the  beam,  the  hydrodynamic  evolution  and  atomic  physics  of  the  plasma, 
as  well  as  the  transport  of  radiation,  must  be  calculated  self- 

consistently  . 

Discussion  of  the  theoretical  model  can  be  separated  for  convenience 
as  fellows:  (a)  hydrodynamics  and  thermal  conduction,  (b)  ionization  and 

atomic  physics,  (c)  radiation  emission  and  transport,  and  (d)  beam  dynamics 
and  energy  deposition. 
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Qualitative  comparisons  of  the  ion  beam-slab  interaction  with  a 
1  7 

treatment  of  a  laser  beam-slab  interaction  shows  that  the  laser  produces 
higher  temperatures  and  lower  densities  on  the  front  (beam)  side  of  the 
plasma,  while  the  proton  beam  is  more  efficient  at  heating  up  the  back  side 
due  to  deeper  penetration  into  the  target. 

Since  these  treatments  of  the  beam-target  interaction  are  one- 
dimer.sicnal ,  energy  flow  normal  to  the  direction  of  the  beam  has  been 
neglected;  in  addition,  self-field  effects  have  not  been  considered.  These 
effects  could  alter  significantly  the  range,  the  magnitude,  and  shapes  of 
the  plasma  profiles  given  here.  A  two-dimensional  radiation  hydrodynamics 
code  will  ce  used  to  investigate  some  of  these  effects  in  a  separate  paper. 


116 


early  times.  In  both  cases,  the  dense  part  of  the  plasma  at  later  times  is 
heated  predominantly  by  radiation  transport  rather  than  thermal 
conduction.  The  laser  case  shows  much  lower  densities  and  higher 
temperatures  on  the  front  side  of  the  plasma,  whereas  on  the  backside  the 
laser  case  shows  much  higher  density  and  lower  temperatures .  Also,  by 
comparing  Figs.  1,2  and  12,  it  can  be  seen  that  the  laser  beam  calculation 
produces  much  sharper  density  and  temperature  gradients.  These  results 
substantiate  the  view  that  proton  beams  are  more  efficient  than  laser  beams 
for  heating  up  the  interior  and  backside  of  the  target;  this  is  due  to 
deeper  penetration  into  the  target  by  the  ion  beams. 

V.  SUMMARY 

The  interaction  of  a  one-MeV  proton  beam  witn  an  aluminum  slab  has 
been  simulated  using  a  fully  self-consistent,  one-dimensional, 
hydrodynamic-ionization-radiation  transport  model.  Detailed  calculations 
of  the  bound  electron  stopping  power  were  used  instead  of  analytic  theories 
such  as  the  L3S  model  and  the  Bethe  equation.  The  atomic  populations  as 
functions  of  density  and  temperature  were  obtained  from  a  collisional- 
radiati ve-equi librium  assumption.  The  CRE  model  encompasses  the  Saha- 
Boltzmann  and  the  corona  models  in  that  it  reproduces  their  results  in 
their  regions  of  validity;  however,  it  is  more  general  in  that  it  is 
applicable  in  intermediate  regimes  in  which  these  models  do  not  apply. 
Radiation  transport  was  calculated  using  a  hybrid  scheme  that  combines  a 
multi-  frequency  formalism  for  the  continuum  with  a  probability-of-escape 
method  for  the  lines.  This  scheme  takes  into  account  the  effect  of  the 
radiation  on  the  atomic  populations,  and  reduces  to  the  photon  diffusion 
approximation  in  regimes  where  that  approximation  is  valid. 

Profiles  of  the  plasma  density  and  temperature  at  various  times  in  the 
plasma  evolution  have  been  presented.  Calculated  spectral  emissions  from 
the  front  and  rear  surfaces  at  various  times  have  also  been  presented. 

Comparisons  with  calculations  omitting  radiation  show  that  radiation 
transport  can  significantly  affect  plasma  evolution,  leading  to 
substantially  different  temperature  and  density  profiles. 


(e.g.,  see  Fig.  8).  Rearside  heating  at  this  time  is  due  primarily  to 
convection.  In  the  first  treatment,  the  beam  is  still  being  stopped  near 
the  center  of  the  plasma,  but  radiation  transport  is  keeping  the  rear  side 
at  a  high  temperature  even  after  the  beam  turns  off  at  10  nsec. 

Fig.  10  shows  the  time  behavior  of  the  density  and  temperature  of  the 
front  side  of  the  plasma.  Here  the  effects  of  radiation  transport  are  more 
pronounced;  the  radiationless  case  produces  much  higher  temperatures  and 
much  lower  densities.  The  marked  drop  in  temperature  after  10  nsec  in  the 
radiation  transport  case  is  due  to  radiative  cooling. 

Fig.  11  gives  the  energy  history  of  the  plasma  along  with  the  energy 
partitioning  for  the  radiationless  case.  The  plasma  energy  is  equal  to  the 
deposited  energy  at  all  times.  The  thermal  energy  is  much  higher  than  in 
Fig.  6  for  the  radiation-hydro  case;  this  indicates  that  radiation  cooling 
occurs  at  the  expense  of  the  plasma  thermal  energy.  This  is  also  indicated 
by  the  much  higher  temperatures  in  the  radiationless  curves  in  Figs.  7-8. 
The  potential  energies  are  not  very  different  in  the  two  calculations.  The 
radiationless  case  has  higher  kinetic  energy  at  later  times;  this  accounts 
for  the  greater  widths  of  the  plasma  in  this  case. 

These  comparisons  of  the  radiationless  and  hydro-radiation  transport 
cases  clearly  establish  the  importance  of  radiation  and  radiation  transport 
effects  in  hot  plasmas  for  the  conditions  considered  here.  Further 
calculations  are  underway  in  order  to  determine  the  parameter  space  over 
which  this  situation  prevails. 


IV.  COMPARISON  WITH  LASER-TARGET  INTERACTION 

It  is  of  interest  to  make  some  qualitative  comparisons  with  a 
simulation  of  a  la3er-aluminum  target  interaction  at  roughly  the  same 
incident  energy  intensity.  Duston,  et .  al .  '  studied  the  interaction  of  a 
1.06  micron  laser  beam  at  an  intensity  of  10  J  W/cm  and  a  long  pulse 
length  (3  nsec  full  width  at  half  maximum)  with  an  8  pm  thick  aluminum 
foil.  Fig.  12  is  taken  from  that  study. 

In  the  laser  case,  after  the  initial  blow  off  plasma  is  formed,  the 
beam  does  not  penetrate  into  the  target  beyond  the  critical  density 
surface,  whereas  the  ion  beam  penetrates  deep  into  the  target,  at  least  at 
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Figure  6  shows  the  energy  history  of  the  plasma  and  the  partitioning 
of  that  energy.  At  early  times,  the  radiated  power  is  an  insignificant 
fraction  of  the  beam  energy,  but  at  late  times,  it  becomes  a  significant 
portion  of  the  energy  budget.  At  about  12  nsec,  it  surpasses  the  total 
energy  (sum  of  kinetic,  potential,  and  thermal  energies)  in  the  plasma.  At 
this  point,  more  than  half  the  energy  deposited  by  the  beam  has  been  lost 
due  to  radiation1  3. 

To  investigate  the  importance  of  radiation  transport,  the  simulation 
was  repeated,  but  with  radiation  omitted.  The  CRE  ionization  dynamics 
treatment  was  retained,  but  the  radiation  was  not  included  in  the  source 
and  sink  terms  in  the  hydrodynamic  equations,  and  optical  pumping  of  the 
atomic  levels  was  excluded.  Fig.  7  shows  temperature  and  density  profiles 
at  3  nsec  with  and  without  radiation.  Including  radiation  effects  produces 
a  more  uniformly- dense  plasma  in  this  case.  The  densities  are  not  very 
different,  except  that  the  profile  without  radiation  shows  a  dip  in  the 
center  and  a  much  less  dense  blow  off  region  in  the  front.  The 
radiationless  temperature  profile  also  attains  a  higher  maximum  value. 

Figure  8  shows  similar  comparisons  at  10  nsec.  Both  density  profiles 
exhibit  gradients  near  the  rear  surface,  but  the  radiationless  density 
gradient  is  much  steeper,  and  has  a  lower  density  blow  off  region  in  the 
front.  Because  there  is  no  energy  loss  due  to  radiation,  much  higher 
temperatures  result. 

Figure  9  shows  the  temperature  and  density  of  the  rearside  of  the 
plasma  as  a  function  of  time  when  radiation  is  included  and  when  it  is 
omitted.  The  curves  are  very  similar  until  about  4  nsec;  after  this  time 
the  rear  side  i3  denser  and  cooler  when  radiation  is  omitted.  Since,  at 
these  times,  the  beam  is  stopped  in  the  plasma  well  before  it  reaches  the 
back  side,  the  large  temperature  difference  between  4  nsec  and  10  nsec 
indicates  that  radiation  transport  plays  a  significant  role  in  heating  the 
rear  side.  Thi3  agrees  with  a  conclusion  reached  in  a  study  of  the 

interaction  of  a  lO1^  W/cm^  laser  pulse  at  1.06  pm  wavelength  on  a  thin 
8.0  pm  A 2,  foil  [Duston,  et.al1^]. 

The  density  and  temperature  rise  around  9  nsec  in  the  radiationless 
case  is  due  both  to  a  shock  reaching  the  rear  surface  and  to  deeper 
penetration  by  the  beam.  This  deeper  penetration  results  from  generally 
lower  plasma  densities  at  later  times  than  when  radiation  is  included 
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The  15  urn  thickness  is  chosen  because  it  corresponds  approximately  to 
the  range  of  1  MeV  protons  in  cold  solid  density  Al^.  Hence,  the  beam  is 
initially  stopped  within  the  target  and  deposits  most  of  its  energy  near 
the  back  surface  where  the  Bragg  peak  occurs  in  the  stopping  cross 
sections . 

A  point  to  be  noted  also  in  this  regard  is  that,  for  consistency,  the 
average  charge  obtained  from  solving  the  CRE  equation  is  used  throughout 
the  calculation.  Continuum  lowering  is  not  taken  into  account.  This  means 
that,  near  solid  density,  the  average  charge  Z,  and  hence  the  electron 
density,  is  probably  low.  The  free  electron  contribution  to  stopping  power 
is  somewhat  underestimated  at  early  times  in  the  plasma  evolution. 
However,  on  a  time  scale  of  the  order  of  nanoseconds,  the  plasma  density  is 
two  orders  of  magnitude  below  solid  density  p  ,  and  the  CRE  Z  is  valid. 

Figures  1  and  2  show  density  and  temperature  profiles  at  different 
times  as  the  plasma  evolves.  As  may  be  noted  from  these  figures,  the 

plasma  spreads  and  becomes  more  uniform  in  density  and  temperature  as  time 
increases . 

Figures  3  through  5  show  calculated  spectra  emitted  from  the  front  and 
rear  sides  of  the  plasma.  In  Fig.  3.  at  0.5  nsec,  the  rear  side  at 

temperatures  near  HO  eV  is  about  twice  as  hot  as  the  front  surface  due  to 

peaking  of  the  stopping  power  near  the  rear  surface.  This  results  in  a 
more  intense  rearside  spectrum.  As  may  be  seen  from  Fig.  2,  the 

temperature  situation  has  reversed  at  later  times,  and  at  4.31  nsec  (Fig. 
4),  the  front  side  is  emitting  a  much  more  intense  spectrum  than  the  rear 
3ide.  The  peak  intensity  on  the  front  side  has  risen  two  orders  of 
magnitude  and  shifted  to  higher  photon  energies.  The  rear  side  intensity 
has  also  increased,  but  less  dramatically.  Some  K  shell  emission  in  the 
region  around  one  to  two  keV  can  be  noted. 

Front  and  rearside  spectra  at  about  10  nsec  (when  the  beam  is  shut 
off)  are  shown  in  Fig.  5.  The  front  side  spectrum  is  very  similar  to  the 
4.31  nsec  spectrum,  but  the  rear  side  intensity  has  increased  as  the  plasma 
continues  heating  and  becomes  more  uniform  [see  Figs.  1  and  2].  K  shell 
emission  i3  very  prominent  here.  At  later  times  when  the  plasma  has  cooled 
considerably  and  become  very  uniform,  the  spectra  from  both  sides  are  very 
similar  and  have  been  reduced  in  intensity  from  the  10  nsec  levels. 
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where  mp  is  the  proton  mass ,  and 
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mB  m 


mB  ♦  m. 
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4.  Enhanced  Stopping  Power  From  Collective  Effects 

The  stopping  power  of  a  dense  ion  beam  can  differ  from  the  sum  of 
single  particle  stopping  powers.  One  way  in  which  this  can  occur  is 
through  phase  mixing  of  the  polarization  wakes  produced  in  the  ambient 
medium  by  the  beam  particles.  The  magnitude  of  this  particular  effect  is 
similar  to  that  given  by  McCorkle  and  Iafrate^2  for  electron  beams, 

S  =  S  (1  +  2ir  n  a3/3)  (28) 

bo  B 

where  Sb  is  the  enhanced  stopping  power,  SQ  is  the  single  particle  stopping 
power,  n3  is  the  beam  number  density  and  a  is  the  ratio  of  the  beam 
velocity  and  local  plasma  frequency 


a 


W 


(29) 
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For  the  case  considered  below,  a  1  MeV  proton  beam  with  a  flux  n0  v0  =  10 
-  ?  —  1 

cm  sec  on  a  solid  density  aluminum  target,  the  correction  to  the 
stopping  power  is  quite  3mall . 


III.  RESULTS 

The  interaction  of  a  monoenergetic  1  MeV  proton  beam  with  a  planar 
aluminum  slab  of  15  ym  thickness  is  investigated.  The  beam  is  assumed  to 
consist  of  a  square-shaped  pulse  of  10  nanosecond  duration  with  a  flux  of 
102^  proton3/cm2-sec.  The  beam  intensity  on  target  is  I.SxlO1^  w/cm2, 
which  is  comparable  to  the  intensities  available  with  current  pulsed-power 
generators . 


Ill 


The  first  term  in  S^(V)  is  the  short  range  ion-electron  binary- 
encounter  scattering  term  and  is  taken  from  the  work  of  Campbell,  7  who 
adapted  it  from  Brueckner  and  Brysk.30  The  second  term  is  the  polarization 
term  and  is  taken  from  Pines  and  Bohm.^  For  distances  larger  than  D,  the 
plasma  acts  as  a  continuous  medium,  and  distant  collisions  cause  loss  of 
energy  by  the  excitation  of  plasma  oscillations,  which  appear  as  an 
oscillating  wake  oehind  the  projectile. 


3.  Interaction  with  Plasma  Ions 


For  stopping  power  due  to  the  plasma  ions,  an  expression  taken  from 
Campbell2^  and  Mehlhorn1^  is  used; 


where 


(dE)_  =  pOO  g  (v) 

dx;i  m.  1 
1 


(21  ) 


w  her  e 


4  2-2 

4ire  Z  Z 

S.(V)  =  - P -  In  A  FU) 

mV 


(22) 


52  =  miV2/2kT.  , 


2 

F(0  =  erf (O  -  4-  O  ♦  =*■)  5®  * 
/ir  mB 


(23) 

(24) 


In  A  -  (1/2)  Zn  (1  +  b  2  /  b  2  ) 

max  min 


(25) 


bmax  is  3et  e(lual  t0  D>  and  bmin  is  8iven  by 


Z  Z  e2 

bmin  =  MAX  ^  _  "2  ; 

umpV  P 


(26) 
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The  bound  electron  stopping  power  is  calculated  from 


(jf)b  -  (pOO/V  sb(v)  (15) 

where  Sb(v)  is  the  stopping  cross  section  in  erg  cm^  and  is  obtained  from 
Eq.  (1H). 

2.  Interaction  With  Free  Electrons 

For  heated  target  materials,  the  atoms  become  ionized,  and  stopping 
due  to  free  electrons  must  be  considered.  The  free  electron  stopping  power 
is  calculated  from 


p  (x) 


Sf(V)  , 


(16) 


with 

4  2  - 

2  ire  Z*  Z  2  . 

Sf(V)  =  - ? -  FU)  {  in  (l  +  +  in  (1  ♦  40  },  (17) 

mV  b 

where  Z  is  the  average  charge  of  the  target  ions,  D  is  the  Debye  shielding 
length,  and 


5  =  (mBV2/2kT)1/2  , 

where  kT  is  the  electron  temperature,  and 

_£2 

F(0  -  erf  (5)  -  2  5  e  4  /  /ir  . 


(18) 


(19) 


The  quantity  b  is  the  minimum  impact  parameter  for  electron-ion 
scattering  and  is  given  by 


b  =  MAX  { 


(20) 


i.e.,  the  maximum  of  either  the  classical  or  quantum-mechanical  impact 
parameter  defined  by  the  uncertainty  principle. 


By  use  of  least-squares  methods,  these  MPP  data  were  fitted  by  a 
function  of  the  form 

Fq(y)  =  exp  (Aq  +  3^y  +  Cqy2  +  Dqy3)  ,  (14) 

where  y  =  In  (E/EQ),  and  EQ  =  0.1  MeV.  F_(y)  gives  stopping  cross  sections 

-ic  o 

in  units  of  10  J  eV  cm  /atom.  Table  1  gives  the  fitting  parameters  Aq , 
Bq,  Cq  and  Dq  for  each  value  of  q.  All  the  fits  were  within  7%,  except  for 
q=3,  where  the  maximum  deviation  was  9%.  Thus,  Eq .  (14)  yielded  reasonably 
good  fits  to  the  calculated  cross  sections. 

In  the  treatment  of  the  beam-target  interaction,  the  average  ionic 
charge  Z  in  a  given  cell  is  used  to  obtain  stopping  power.  For  non-integer 
values  of  2,  Eq .  (14)  is  used  for  q  and  q+1  ,  where  q  is  the  integer  part 
of  Z,  and  the  stopping  cross  section  for  Z  is  obtained  by  linear 
interpolation  between  Fq(y)  and  Fq+1(y). 

For  q=0,  q=*1  and  q»2,  the  MPP  data,  and  hence  the  fits,  are 

monotoni cally  decreasing  for  energies  >_  0.1  MeV.  For  these  cases,  a  local- 
oscillator  model  (LOM)2"’2®,  is  used  to  calculate  bound  electron  stopping 
power  for  energies  below  0.1  MeV.  These  LOM  cross  sections  are  adjusted  to 
match  the  MPP  data  at  0.1  MeV,  and  the  cross  sections  at  lower  energy  are 
scaled  accordingly.  The  MPP  data  for  q  »  8,  9,  and  10  are  also  monotically 
decreasing  for  E  >_  0.1  MeV,  but  at  these  charge  levels,  the  stopping  power 
is  dominated  by  the  free-electron  contri bution;  hence,  no  adjustments  to 
the  bound  electron  stopping  power  at  energies  less  than  0.1  MeV  are  made. 
The  fits  are  just  extended  below  0.1  MeV. 

For  the  other  q  values  between  3  and  11,  the  MPP  data  peak  in  the  0.1 
MeV  -  100  MeV  range;  hence  the  fits  are  extended  below  0.1  MeV  to  obtain 
low  energy  cross  section  data.  For  q  >  11,  the  bound  electron  stopping 
power  is  set  to  3ome  low  value,  since  free  electron  stopping  dominates. 

From  an  energetics  point  of  view,  errors  in  the  low  energy  regime  are 
not  as  critical  as  at  higher  energies.  For  example,  a  proton  with  an 
initial  energy  of  one  MeV  has  lost  90%  of  its  energy  and  is  near  the  end  of 
its  range  by  the  time  it  degrades  to  a  0.1  MeV  proton.  Therefore,  errors 
in  stopping  power  below  0.1  MeV  wili  not  affect  energy  deposition  profi.es 
very  much. 


108 


d.  Beam  Dynamics  and  Energy  Deposition 


It  is  assumed  that,  on  the  timescale  of  the  slab  response  (several 
nanoseconds),  the  beam  maintains  uniform  flow,  that  is 


=  m 


(10) 


where  p0  and  Ug  are  the  local  mass  density  and  axial  velocity  of  the 
beam.  The  mass  flux  m  is  assumed  constant  up  to  the  point  the  beam  is 
stopped.  J  is  the  beam  current  per  unit  area;  mg  and  2g  are  the  mass  and 
charge  of  a  beam  ion. 


For  a  relativistic  beam  with  negligible  thermal  spread,  the  kinetic 
energy  of  a  beam  ion  is  related  to  the  local  velocity  by 

2 


mgc 


JB 


Al-Ug/c2) 


-  mBc 


(11) 


The  energy  is  also  related  to  the  stopping  power  by 


E 


B 


(x  ) 


E 


o 


(-) 

'•dx'T 


ax 


(12) 


where  EQ  is  the  initial  energy,  x'  is  the  distance  into  the  slab,  and  the 
total  stopping  power  is  given  by  the  sum  of  the  ion-free  electron,  ion- 
bound  electron  and  ion-plasma  ion  interactions, 


r dE i  _  fdE'i  +  f^.1  +  f— ) 
^dx^T  =  ^dx^e  Mx'b  +  ^dx'i 


(13) 


1.  Interaction  With  Bound  Electrons 

Recently,  McGuire,  Peek,  and  Pitchford1^  (MPP)  published  detailed 
calculations  of  stopping  power  cross  sections  of  protons  by  Ai  ions.  By 
using  the  generali zed-oscillator-strength  formulation  of  the  Born 
approximation,  ‘they  calculated  tables  of  stopping  cross  sections  for 
protons  interacting  with  At  ions  of  charge  q,  where  0  £  q  £  11,  in  the 
energy  range  0.1  MeV  to  100  MeV. 


P  u 

The  multif requency  model 41  solves  the  equation  of  radiation  transport 
at  a  large  number  of  discrete  frequencies,  providing  resolution  of  emission 
lines,  recombination  edges  and  absorption  edges.  It  provides  accurate 
radiation  transport  at  high  density,  and  gives  spectral  details  such  as 
self-absorption  features.  However,  a  large  number  of  frequencies  is 
required  to  provide  adequate  resolution,  and  it  tends  to  be  a  relatively 
costly  method. 

The  hybrid  model  uses  the  multifrequency  formalism  to  transport  the 
continuum  (bound-free  and  free-free)  radiation  and  frequency-integrated 
escape  probabilities  to  transport  the  lines.  Continuum  opacities  are 
interpolated  from  the  multifrequency  mesh  and  folded  into  the  line 
transport  calculations.  The  line  opacities  are  assumed  to  have  negligible 
impact  on  the  continuum  energetics. 

In  all  of  the  models  outlined  above,  local  ionization  state- dependent 
inner-shell  opacities  are  included,  since  these  processes  are  very 
important  in  the  cool,  dense  plasma  regions.  Inner-shell  photoionization 
cross  sections  for  the  neutral  element  are  taken  from  the  fits  by  Biggs  and 
Lighthill2^.  For  ions,  inner  shell  photoionization  cross  sections  are 
estimated  using  the  neutral  cross  section  data,  but  with  the  absorption 
edges  shifted;  the  positions  of  the  ioni zat ion- dependent  absorption  edges 
are  based  on  the  Hartree-Fock  calculations  of  dementi  and  Roetti.2^ 

The  local  rate  of  energy  change  in  zone  j ,  due  to  radiation  transport 
is  given  by 

;j  -  - 1  <fpj  - 1  cpkj  '«>  (<» 

where  Fpi<  is  the  rate  of  energy  loss  in  zone  k  due  to  a  discrete  radiative 
process  (or  frequency  group)  P,  and  Cpkj  is  the  radiative  coupling  of  zone 
k  to  zone  j  for  that  process.  The  couplings  are  functions  of  opacity, 
integrated  over  process  and  photon  path.  In  the  probabilistic  model,  a 
matrix  of  couplings  must  be  computed  for  each  bound-bound,  bound-free  and 
free-free  process;  for  the  multifrequency  model,  they  must  be  computed  for 
each  discrete  frequency.  In  this  way,  the  net  cooling  and  heating  by 
radiation  emission  and  absorption  among  the  various  zones  of  the  plasma  is 
accounted  for  accurately. 
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The  ionization  and  excitation  energy  can  also  be  calculated  by 


Z  x- 
.1  i 
1 


I  ’ 


(8) 


where  is  tfie  energy  of  level  i,  measured  from  the  ground  state  of  the 
neutral  atom. 


c.  Radiation  Emission  and  Transport 

Radiation  emission  from  a  plasma  and  its  opacity  are  dependent  on  the 

local  atomic  level  population  densities.  Except  for  optically  thin 

plasmas,  however,  the  level  populations  depend  on  the  radiation  field, 

since  optical  pumping  via  photoionization  and  photoexcitation  can  produce 

significant  population  redistribution.  Thus,  the  ionization  and  radiation 

transport  processes  are  strongly  coupled  and  must  be  solved  self¬ 
op 

consistently.  In  this  model,  an  iterative  procedure  is  used,  where  level 
populations  are  calculated  using  the  radiation  field  from  the  previous 
iteration,  then  using  these  populations  to  calculate  a  new  radiation  field 
and  recalculating  populations  until  convergence  is  reached. 

Three  distinct  radiation  transport  schemes  have  been  developed,  and 
can  be  used  interchangeably  in  the  code:  A  probabilistic  scheme,  a 

multifrequency  scheme  and  a  hybrid  method.  Descriptions  of  these  transport 
methods  have  appeared  elsewhere;  in  this  paper,  their  general  properties 
will  be  outlined. 

The  probabilistic  model forms  local  angle  and  frequency  averaged 
escape  probabilites  for  each  emission  line  and  for  each  bound-free 
process.  Free-free  radiation  is  treated  with  a  multifrequency  formalism. 
The  radiation  transport  and  emission  spectra  are  calculated  from  these 
escape  probabilities.  The  method  is  cost-effective,  can  treat 

comprehensive  atomic  models  and  provides  good  overall  energetics,  but 
cannot  calculate  certain  spectral  details  and  breaks  down  at  very  high 
densities . 


thermal  equilibration  times  are  short  compared  with  the  time  scales  of  the 
energy  input  and  changes  in  the  hydrodynamic  variables. 

The  local  rate  of  change  of  energy  due  to  radiation  transport  e  , 

i  du 

and  that  due  to  the  beam  deposition  edep»  will  be  discussed  in  Sections  C 
and  D . 


o.  Ionization  and  Atomic  Physics 


The  ionic  populations  in  the  plasma  may  be  characteri zed  by  a  set  of 
atomic  rate  equations  of  the  form 


df . 
1 

dt 


Z .  W..  f. 
j  Ji  J 


z.  w. .  f. 

j  1J  1 


v6) 


where  f^  is  the  fractional  population  of  atomic  level  i,  and  is  the  net 
reaction  rate  describing  the  transition  from  initial  state  j  to  final  state 
i.  An  equation  of  this  type  is  constructed  for  each  of  the  atomic  levels 
included  in  the  model. 


For  sufficiently  dense  plasmas  of  the  sort  we  expect  to  model  with  the 
ion-beam  code,  the  effective  populating  and  depopulating  rates  are 
generally  fast  compared  with  the  hydrodynamic  response.  An  equilibrium 
assumption  can  be  justified,  which  involves  dropping  the  explicit  time 
dependence  in  equation  (6).  The  plasma  is  then  said  to  be  in  collisional- 
radiative  equilibrium  (CRE),  whereby  the  plasma  ionization  state  responds 
instantaneously  to  changes  in  hydrodynamic  quantities. 


The  rate  coefficients  that  are  used  to  calculate  the  populating  and 
depopulating  rates,  wji»  are  calculated  using  various  scattering 
techniques.  The  processes  included  in  this  calculation  and  the  methods 
used  in  calculating  the  corresponding  rate  coefficients  are  summarized 
elsewhere  ' . 

Once  the  set  of  rate  equations  (including  the  radiation  transport)  has 
been  solved  for  the  level  populations,  f ^ ,  the  electron  density  can  be 
calculated , 


N 


e 


E  zi  fi  NI 


(7) 


where  z ^  is  the  ionic  charge  of  level  i  and  is  the  total  ion  density. 
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a.  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  mass,  momentum,  and  total  energy 

are  transported  in  one  dimension  using  a  numerical  scheme  with  a  sliding- 

20 

zone  version  of  flux-corrected  transport.  A  special  gridding  algorithm 
is  used  which  moves  zones  in  a  Lagrangian  fashion  and  adjusts  the  mesh  in 
order  to  resolve  steep  gradients  in  the  flow.  The  hydrodynamic  equations 
solved  are 


Du  3p  3  ,  >  - 

Dt  a  at  +  37  (up)  =  °* 


0(pu) 

Dt 

De 


d P 
dx  ’ 


T  3  •  • 

— -  =s  -  —  (up)  +  e  +  e 
Dt  3x  K  J  Erad  ^dep 


*  h  1  "nht  1  • 


(1) 

(2) 


(3) 


where  o  is  mass  density,  u  is  velocity,  F  is  pressure,  eT  is  total  energy 
density,  erad  is  the  rate  of  energy  loss  or  gain  due  to  radiation, 
is  the  rate  of  energy  gain  due  to  the  beam  deposition  process,  n  is  the 
thermal  conducti vi  ty ,  and  N  is  the  ion  density.  The  thermal  conduction  is 
calculated  implicitly,  using  an  iterative  Crank-Nicholson  scheme. 

Since  densities  did  not  much  exceed  solid  density  in  this  study,  a 
simple  equation  of  state  was  assumed  with 

P  =  |  (eT  -  1  pu2  -  ej)  ,  (4) 

where  is  the  ion  potential  energy  due  to  ionization  and  excitation.  A 
single  temperature  model  was  employed, 

kT  =  - - -  ,  (5) 

(p/mj)  (1+Z) 

where  trij  is  ion  mass,  and  T  is  temperature.  The  ionization  energy,  e^.,  and 
effective  charge,  Z  are  calculated  from  the  ionization-radiation  equations 
which  are  explained  below.  A  single  temperature  assumption  is  valid  in  the 
solid  dense  material  as  well  as  in  the  beam  deposition  region,  since 


r* . 
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Table  1  —  Least-square-fitting  parameters  to  the  MPP  stopping  powers* 


A 

q 

B 

q 

C 

q 

D 

q 

3.135 

-4.415 

(-1) 

-2.592 

(-2) 

-1.031 

(-3) 

2.773 

-3.534 

(-D 

-3.719 

(-2) 

-5.234 

(-4) 

2.159 

-8.853 

(-2) 

-9.409 

(-2) 

3.625 

(-3) 

1.012 

6.257 

(-1) 

-2.644 

(-1) 

1.651 

(-2) 

8.459 

(-D 

5.977 

(-1) 

-2.480 

(-1) 

1.482 

(-2) 

6.575 

(-1) 

5.506 

(-1) 

-2.331 

(-1) 

1.374 

(-2) 

5.423 

(-1) 

4.501 

(-1) 

-2.085 

(-1) 

1.212 

(-2) 

4.447 

(-1) 

2.828 

(-D 

-1.564 

(-1) 

7.782 

(-3) 

3.716 

(-1) 

8.571 

(-2) 

-1.072 

(-1) 

4.287 

(-3) 

3.148 

(-1) 

-1.870 

(-1) 

-3.894 

(-2) 

-4.858 

(-4) 

-4.684 

(-1) 

-2.246 

(-1) 

6.746 

(-3) 

-5.082 

(-3) 

-7.250 

2.867 

-5.093 

(-1) 

2.388 

(-2) 

Numbers  in  parentheses  indicates  powers  of  ten 
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Fig.  3  —  Front  and  rear  side  spectra  at  0.5  nsec 
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Fig.  6  —  Beam-target  energy  partitioning  as  a  function  of  time,  including  beam 
deposition  (dep),  total  plasma  energy  (sum),  plasma  kinetic  energy  (kin), 
plasma  thermal  energy  (int),  ionization  energy  (pot)  and  radiated  energy  (rad). 


DENSITY  (gm/cm  )  TEMPERATURE  (e> 


r(cm) 

Fig.  7  —  Comparison  of  temperature  and  density  profiles  at  3  nsec 
with  and  without  effects  of  radiation 
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Fig.  8  —  Comparison  of  temperature  and  density  profiles  at  10  nsec 
with  and  without  effects  of  radiation 
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Fig.  9  —  Temperature  and  density  history  of  the  rear  surface  of  the  plasma 
with  and  without  radiation  transport.  T  is  temperature  and  p  is  density. 
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Fig.  11  —  Beam-target  energy  partitioning  when  radiation  effects  are  omitted  (same  labels  as  Fig.  6). 
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Fig.  12  —  Temperature,  density,  and  radiative  cooling  profiles  at  about  5  nsec 
in  the  laser  heated  foil  study  (Ref.  31). 
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